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Abstract 

This  report,  prepared  by  members  of  the  Working  Group  on  Thermal  Regimes 
of  the  Section  on  Ice  Research  and  Engineering  of  the  International  Association 
for  Hydraulic  Research,  is  a  comprehensive  overview  of  frazil  ice.  Starting  from 
the  initial  nucleation  of  single  frazil  ice  crystals  to  the  formation  of  ice  covers  that 
may  be  many  kilometers  long,  the  physics  and  evolution  of  frazil  ice  in  natural 
water  bodies  are  described.  Laboratory  experiments  conducted  over  the  last  30 
years  on  frazil  ice  dynamics  and  other  aspects  of  frazil  are  described  and 
classified.  A  physically  based,  quantitative  model  that  describes  the  dynamic 
evolution  of  the  crystal  size  distribution  function  is  presented.  In  addition,  the 
development  of  numerical  models  of  frazil  ice  in  oceans  and  rivers  is  discussed 
and  their  results  described .  This  report  serves  as  a  review  of  the  state  of  the  art 
of  the  present  understanding  of  frazil,  and  the  extensive  references  are  a 
comprehensive  resource. 
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CHAPTER  1 
Physics  of  Frazil  Ice 

THOMAS  FOREST 


NOMENCLATURE 


NUCLEATION  OF  FRAZIL  CRYSTALS 
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molar  concentration  dissolved  in 
liquid 

latent  heat  of  fusion 

universal  gas  constant 

critical  radius 

interface  temperature 

liquid  temperature 

freezing  point  of  liquid 

freezing  point  of  pure  liquid 
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xm  1=x> 

interface  supercooling 
T  -T 

specific  volume  of  solid 
growth  velocity 

solid/liquid  interface  roughness 

surface  free  energy 

fraction  of  binding  energy  associated 

with  crystal  face 

constants 


INTRODUCTION 

The  formation  and  evolution  of  "frazil"  ice  in  tur¬ 
bulent  water  that  undergoes  a  small  degree  of  su¬ 
percooling  is  a  problem  that  involves  several  funda¬ 
mental  aspects.  The  important  questions  regarding 
the  basic  physics  of  the  process  are: 

1.  How  do  frazil  ice  crystals  originate? 

2.  What  shape  do  frazil  crystals  adopt  and 
why? 

3.  How  is  the  size  of  frazil  crystals  related  to 
the  thermodynamic  concept  of  critical  ra¬ 
dius? 

In  this  chapter,  we  present  some  current  ideas  about 
these  questions. 


Nucleation  is  a  general  term  referring  to  the  for¬ 
mation  of  a  new,  embryonic  phase  from  a  parent 
phase.  In  this  case,  ice  crystallizes  from  water,  which 
is  the  parent  phase.  The  fundamental  thermody¬ 
namic  requirement  for  nucleation  to  take  place  is  a 
certain  level  of  supercooling  in  the  parent  phase,  i.e., 
the  temperature  of  the  water  must  be  below  the 
equilibrium  freezing  point  of  water;  if  the  water  is 
pure  (no  dissolved  impurities)  the  freezing  point  is, 
by  definition,  0°C  at  atmospheric  pressure.  If  the 
water  is  indeed  pure  and  contains  neither  dissolved 
impurities  nor  any  undissolved  particles  such  as 
dirt,  bacteria,  etc.,  then  ice  nucleation  will  not  occur 
until  an  impressively  large  supercooling  is  attained. 
This  limit  is  called  the  homogeneous  nucleation  su¬ 
percooling,  which  for  water  is  in  the  neighborhood 
of  40°C,  i.e.,  the  water  must  be  cooled  to  -40°C  be¬ 
fore  ice  will  nucleate  (Mossop  1955). 

Natural  bodies  of  water  are  not  pure;  they  con¬ 
tain  both  dissolved  and  undissolved  materials.  Un¬ 
der  these  conditions,  the  undissolved  particles  can 
act  as  a  catalyst,  whereby  ice  nuclei  can  form  on 
these  particles  at  supercoolings  that  are  much 
smaller  than  the  homogeneous  limit.  This  type  of 
nucleation  is  called  heterogeneous  nucleation.  The 
ability  of  these  particles  to  enhance  nucleation  is  re¬ 
lated  to  the  geometry  and  chemistry  of  their  sur¬ 
faces.  Perhaps  the  best  known  ice  nucleating  cata¬ 
lyst  is  silver  iodide,  which  is  an  inorganic  solid  and 
has  a  heterogeneous  nucleation  supercooling  of  3  to 
4°C  (Fletcher  1968).  Organic  particles  can  also  be  ef¬ 
ficient  ice  nucleators;  certain  bacteria  can  nucleate 
ice  at  supercoolings  as  low  as  1°C  (Lindow  et  al. 
1978).  Unfortunately,  natural  bodies  of  water  sel¬ 
dom,  if  ever,  are  supercooled  by  1°C.  Measurements 
in  streams  just  prior  to  the  onset  of  frazil  ice  growth 
indicate  a  much  smaller  supercooling,  on  the  order 
of  0.1°C  or  less  (Carstens  1966,  Osterkamp  et  al. 


1973).  At  this  level  of  supercooling,  there  is  no 
known  catalyst  that  has  such  a  low  heterogeneous 
nucleation  limit.  This  has  led  to  the  concept  of  sec¬ 
ondary  nucleation  or  collision  breeding  (Strickland- 
Constable  1972,  Botsaris  1976)  as  a  mechanism  for 
the  creation  of  frazil  ice  crystals. 

Secondary  nucleation  refers  to  the  process  where 
tiny  ice  fragments  (on  the  order  of  a  few  microme¬ 
ters  in  size)  are  shed  from  an  existing  large  ice  crys¬ 
tal  when  it  collides  with  a  solid  surface  in  the  water 
or  when  it  collides  with  another  large  ice  crystal.  Ex¬ 
periments,  using  both  ice  and  salt  crystals  (Garabed- 
ian  and  Strickland-Constable  1974,  Garside  and  Lar¬ 
sen  1978),  have  indicated  that  large  numbers  of 
small  fragments  are  generated  by  low-energy  colli¬ 
sions  of  crystals  at  low  supercoolings.  These  tiny 
fragments  then  act  as  growth  centers  for  new  frazil 
ice  crystals.  In  this  way,  a  few  seed  crystals  can  very 
quicldy  generate  first,  second,  third,  etc.,  generation 
ice  nuclei  that  grow  into  frazil  ice  crystals,  as  long  as 
supercooling  persists  in  the  water.  It  is  easy  to  see 
that  secondary  nucleation  can  quickly  generate  the 
number  of  frazil  ice  crystals  that  are  often  observed 
in  laboratory  and  field  studies  (frazil  ice  concentra¬ 
tions  are  typically  on  the  order  of  106  crystals /m3 
[Schaefer  1950]).  The  two  conditions  necessary  for 
secondary  nucleation  are  1)  turbulent  motion  of  the 
water,  and  2)  the  presence  of  seed  ice  crystals. 

Sufficient  turbulent  motion  is  generated  in 
streams  and  rivers  where  the  current  is  greater  than 
approximately  1.0  m/s  (Carstens  1970);  in  lakes,  sur¬ 
face  turbulence  can  be  created  by  wind  and  waves. 

With  turbulence  and  a  small  degree  of  supercooling, 
frazil  ice  will  be  produced  if  some  seed  ice  crys¬ 
tals  are  introduced  into  the  water.  These  seed 
crystals  can  come  from  a  number  of  different 
sources.  Since  frazil  ice  generally  forms  when 
ambient  air  temperatures  are  lower  than  approxi-  -  - 

mately  -8°C,  vapor  evaporating  from  the  water 
surface  will  encounter  cold  air  and  can  sublimate 
into  ice  crystals.  These  ice  crystals  will  then  fall 
down  to  the  water  surface  where  they  are  en-  0.741 
trained  by  the  turbulent  motion  of  the  water.  nm 
Osterkamp  et  al.  (1974)  have  observed  this  seed¬ 
ing  mechanism  in  the  field. 

Small  water  droplets  are  also  generated  at  the 
surface  of  turbulent  water  by  breaking  waves 
and  bubbles  bursting  at  the  surface.  These  drop¬ 
lets  are  ejected  into  the  cold  air  where  they  can 
readily  freeze  and  subsequently  fall  back  into  the 
water  to  act  as  seed  crystals  (Gosink  and 
Osterkamp  1986).  Other  possible  sources  of  seed 
crystals  are  border  ice,  which  generally  forms 
along  the  banks  of  a  river,  cold  dirt  particles  that 


fall  to  the  surface  of  the  water  and  act  as  heteroge¬ 
neous  ice  nucleators  (termed  contact  nucleation 
[Cooper  1974]),  and  snow  or  sleet.  Undoubtedly,  the 
seeding  mechanism  is  some  combination  of  the 
above  sources,  depending  on  the  meteorological 
conditions;  however,  once  seeded,  secondary  nucle¬ 
ation  is  responsible  for  the  prodigious  volume  of 
frazil  ice  that  is  often  generated  in  rivers  at  freezeup. 

FRAZIL  ICE  MORPHOLOGY 

Of  the  many  different  ice  crystal  shapes  that  are 
observed  in  nature,  none  are  as  elegant  in  their  sim¬ 
plicity  as  frazil  ice  crystals.  During  growth,  frazil 
generally  evolves  as  thin,  circular,  disk-shaped  crys¬ 
tals.  Depending  on  the  conditions  under  which 
frazil  crystals  grow,  their  maximum  disk  diameter  is 
on  the  order  of  1  or  2  mm,  while  their  thickness  var¬ 
ies  from  approximately  10  to  100  pm.  This  disk  mor¬ 
phology  directly  results  from  the  strong  anisotropic 
growth  rates  for  ice's  two  crystallographic  direc¬ 
tions. 

When  water  freezes,  water  molecules  arrange 
themselves  into  hexagonal  unit  cells  that  act  as  the 
basic  building  blocks  for  all  ice  crystals  (Fig.  1).  For 
our  purposes,  there  are  two  principal  directions  of 
growth:  along  the  c-axis,  which  is  the  axis  to  the  unit 
cell,  and  the  fl-axis,  which  is  perpendicular  to  the  c- 
axis.  The  ice  growth  rates  in  these  two  directions  are 
controlled  by  two  very  different  mechanisms.  Jack- 
son  (1958)  and  Jackson  et  al.  (1967)  have  presented  a 
general  theory  to  characterize  the  roughness  of  a 


nm 


Figure  1.  Shape  of  unit  cell  for  hexagonal  ice. 
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solid /liquid  interface.  They  define  the  parameter  a 
for  a  crystallographic  plane  as 

a  =  -£^  0) 

RT0 

where  e  =  fraction  of  the  total  binding  energy  of 
an  atom  that  can  be  associated  with  an 
atomic  layer  parallel  to  the  face  under 
consideration 
L  =  latent  heat  of  fusion 
R  =  universal  gas  constant 

T0  =  melting  temperature  of  ice. 

Jackson  concluded  that  a  crystal  face  will  be  smooth 
when  a  >  2  and  will  be  rough  when  a  <  2.  The  kinet¬ 
ics  (i.e.,  the  attachment  of  water  molecules  to  the  ice 
surface)  for  a  smooth  crystal  is  more  likely  to  be  so 
slow  that  it  will  control  the  growth  rate;  in  this  case, 
the  crystal  face  will  have  large  flat  faces.  Conversely, 
for  rough  crystal  faces,  the  attachment  kinetics  is 
rapid  enough  that  the  surface  grows  continuously; 
in  this  case,  the  dissipation  of  latent  heat  from  the 
interface  limits  the  growth  rate.  Fujioka  (1978)  has 
shown  that  values  for  a  for  the  a-  and  c-axis  direc¬ 
tions  are  0.88  and  2.64  respectively.  This  implies  that 
the  growth  in  the  a-axis  is  limited  by  the  diffusion  of 
latent  heat,  while  c-axis  growth  is  limited  by  some 
kinetic  attachment  mechanism. 

For  frazil  ice  crystals,  the  c-axis  grows  by  the  con¬ 
secutive  nucleation  of  thin  monolayers  spreading 
across  the  flat  crystal  face.  For  growth  to  occur,  a 
stable  island  must  be  nucleated  on  the  growing  sur¬ 
face  that  quickly  spreads  across,  covering  the  crystal 
face.  The  general  form  for  c-axis  growth  is 

V  =  \ixe^/ATi  (2) 

where  V  =  growth  velocity 
Hi,  P2  =  constants 

AT,  =  interfacial  supercooling,  i.e.,  the  differ¬ 
ence  between  the  interface  temperature 
Tj  and  the  temperature  in  water  far 
away  from  the  interface  TTO. 

For  the  a-axis,  the  general  form  of  growth  rate  is 

V  =  »3ATi  (3) 

where  |X3  is  a  constant  related  to  a  heat  transfer  coef¬ 
ficient.  In  general,  the  value  of  Tj  is  different  for  each 
face  of  the  crystal.  Experiments  (Hillig  1958, 
Michaels  et  al.  1966)  have  shown  that  at  the  low  lev¬ 
els  of  supercooling  under  which  frazil  crystals 
evolve,  the  growth  rate  in  the  c-axis  direction  is  one 
to  two  orders  of  magnitude  smaller  than  that  in  the 


a-axis  direction.  As  a  result  of  these  two  different 
growth  mechanisms,  frazil  crystals  evolve  as  thin, 
circular  disks  with  thickness-to-diameter  aspect  ra¬ 
tios  in  the  range  of  0.05  to  0.2  (Bukina  1963). 


STABILITY  OF  FRAZIL  CRYSTALS 


Frazil  ice  crystals  most  likely  originate  from  the 
tiny  ice  fragments  that  are  generated  by  crystal- 
crystal  collisions.  Even  though  there  is  a  supply  of 
small  ice  nuclei,  a  question  remains  as  to  whether 
these  nuclei  will  grow  or  decay;  furthermore,  if  nu¬ 
clei  grow,  is  there  any  size  limit  to  their  growth?  In 
this  section,  the  thermodynamics  of  crystal  evolu¬ 
tion  in  supercooled  water  will  be  briefly  reviewed. 

The  thermodynamic  concept  of  a  critical  radius  is 
used  to  decide  whether  nuclei  grow  or  dissolve.  For 
any  phase  transition,  a  critical  radius  can  be  defined 
as  the  size  of  a  nucleus  that  is  in  thermodynamic 
equilibrium  with  the  surrounding  parent  phase.  It 
can  be  shown  that  this  equilibrium  state  is  unstable. 
That  is,  if  the  size  of  a  nucleus  is  less  than  the  critical 
size,  it  will  dissolve  back  into  the  parent  phase;  if  it  is 
larger  than  the  critical  size,  it  will  grow  (Turnbull 
and  Fisher  1949).  Note  that  the  critical  size  is  a  ther¬ 
modynamic  variable  just  as  temperature  or  pressure 
is  a  thermodynamic  variable.  For  crystals  in  a  super¬ 
cooled  melt,  the  critical  size  or  radius  Rc  is  given  by 
the  Gibbs-Thomson  equation  (Dufour  and  Defay 
1963) 


where  y  =  surface  free  energy  of  a  solid/liquid  in¬ 
terface 

vs  =  specific  volume  of  the  solid 
L  =  specific  latent  heat  of  fusion 
T0  =  freezing  point  of  the  pure  liquid  (0°C  for 
water). 

The  supercooling  AT  is  the  difference  between  the 
equilibrium  freezing  point  of  the  liquid  Tm  and  the 
actual  temperature  of  the  liquid  TM.  For  most  natu¬ 
ral  bodies  of  water  that  contain  some  dissolved  sol¬ 
ids,  Tm  will  be  slightly  less  than  T0.  This  freezing 
point  depression  ATm  =  T0-Tm  is  given  as 

ATm=RT02C0  (5) 

where  C0  is  the  molar  concentration  of  material  dis¬ 
solved  in  the  liquid  (kmol  dissolved  material/kmol 
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liquid);  eq  5  assumes  that  CG  «  1.  If  the  dissolved 
material,  e.g.,  salts,  dissociate  in  solution,  the  expres¬ 
sion  for  A Tm  is  more  complicated.  Osterkamp  (1978) 
gives  a  thorough  discussion  of  this  point  for  river 
water  and  concludes  that  A Tm  is  on  the  order  of 
0.005°C.  Equation  4  assumes  that  the  nucleus  is 
spherical;  if  the  nucleus  has  some  other  shape,  such 
as  a  disk,  the  expression  for  the  critical  radius  will 
differ  from  eq  4.  For  example,  if  a  frazil  disk  is  mod¬ 
eled  as  a  "stubby"  cylinder  (Fujioka  and  Serkerka 
1974)  then 


If  the  temperature  of  the  liquid  phase  is  equal  to  the 
equilibrium  freezing  point  Tm,  then  AT  approaches 
zero  and  the  critical  radius  becomes  infinite.  Physi¬ 
cally,  this  implies  that  all  nuclei,  regardless  of  their 
size,  will  start  to  dissolve  into  the  liquid  phase. 
Thus,  frazil  ice  crystals  will  only  be  able  to  survive 
and  grow  if  the  water  temperature  is  less  than  the 
freezing  point  Tm. 

The  traditional  derivation  of  the  Gibbs-Thomson 
equation  assumes  that  the  presence  of  a  crystal 
nucleus  does  not  alter  the  thermodynamic  state  of 
the  parent  phase,  i.e.,  supercooling  or  concentration 
of  dissolved  material.  This  is  true  if  the  nucleus  is 
surrounded  by  an  infinite  volume  of  liquid  phase. 
However,  frazil  ice  crystals  appear  in  large  concen¬ 
trations  (on  the  order  of  106/m3),  which  implies  that 
the  thermodynamic  state  of  the  water  will  be  af¬ 
fected  by  their  presence.  For  a  sufficiently  large  su¬ 
percooling  (on  the  order  of  0.1°C),  a  large  fraction  of 
frazil  nuclei  survive  and  begin  to  grow,  since  the  size 
of  the  crystals  is  greater  than  the  critical  size  (which 
is  on  the  order  of  0.40  pm  for  this  level  of  supercool¬ 
ing).  As  each  crystal  grows,  the  ice  phase  rejects  dis¬ 


solved  material  to  the  remaining  liquid  phase.  The 
concentration  of  dissolved  material  in  the  liquid 
phase  slowly  increases,  resulting  in  an  increase  in 
the  equilibrium  freezing  point  depression,  accord¬ 
ing  to  eq  5.  If  this  change  in  concentration  takes 
place  at  a  fixed  liquid  temperature  T,M  then  the  su¬ 
percooling  AT,  defined  above,  decreases.  According 
to  eq  4,  the  critical  size  therefore  increases,  but  re¬ 
mains  less  than  the  physical  size  of  the  crystals.  At 
some  point,  the  increase  in  concentration  of  the  liq¬ 
uid  phase  is  enough  to  make  the  critical  size  equal  to 
the  physical  size  of  the  crystals.  The  crystals  now 
reach  a  new  equilibrium  state  that  is  stable. 

Such  increases  in  concentration  of  the  liquid 
phase  during  periods  of  frazil  growth  have  been 
measured  in  field  studies  by  Osterkamp  et  al.  (1975). 
Such  a  stability  analysis  was  carried  out  (Forest 
1986)  for  frazil  crystals  with  a  distribution  in  crystal 
thickness.  The  thermodynamic  model  predicted 
that  the  stable  equilibrium  radius  of  the  frazil  crys¬ 
tals  increased  with  an  increase  in  crystal  thickness. 
The  results  showed  that  for  the  observed  range  of 
thickness  (approximately  10  to  80  pm),  the  stable 
equilibrium  radii  varied  from  approximately  50  pm 
to  several  millimeters,  depending  on  the  supercool¬ 
ing.  One  of  that  work's  principal  conclusions  was 
that  thinner  disks  will  quickly  reach  a  limiting  size 
that  can  be  relatively  small  and  then  cease  to  grow. 
For  thicker  disks,  the  limiting  size  is  very  large, 
which  implies  that  these  disks  are  free  to  grow  as 
long  as  there  is  enough  supercooling  to  drive  the 
growth  process.  This  type  of  growth  rate  dispersion 
phenomenon,  where  small  crystals  have  a  much 
smaller  average  growth  rate  than  larger  crystals,  has 
been  observed  previously  (Cise  and  Randolph  1972). 

To  completely  understand  the  evolution  of  frazil 
crystals,  we  must  couple  the  thermodynamic  stabil¬ 
ity  analysis  to  an  accurate  model  for  predicting  the 
growth  rates  of  frazil  crystals. 
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CHAPTER  2 

Experimental  Investigation  of  Frazil  Ice 

KATHLEEN  D.  WHITE  AND  STEVEN  F.  DALY 


NOMENCLATURE 

d  pipe  diameter 

/  friction  factor 

Q  heat  loss  rate 

Rh  hydraulic  radius 

ts  time  to  maximum  supercooling 

u  average  velocity 

u*  friction  velocity 

e  turbulent  energy  dissipation  rate 

k  Von  Karman's  constant 

v  kinematic  viscosity 

INTRODUCTION 

Laboratory  experiments  have  been  a  valuable 
tool  for  developing  an  understanding  of  some  as¬ 
pects  of  the  growth  and  evolution  of  frazil  ice  crys¬ 
tals.  As  described  in  Chapter  3,  frazil  ice  can  evolve 
through  many  forms  in  natural  water  bodies.  In 
general,  however,  the  distances  are  too  large  and  the 
time  scale  too  long  to  replicate  this  evolution  in  any 
laboratory  setting.  As  a  result,  laboratory  experi¬ 
ments  have  concentrated  upon  the  formation  phase, 
which  is  characterized  by  supercooled  water,  rap¬ 
idly  growing  disk-shaped  crystals,  and  the  creation 
of  new  crystals  by  secondary  nucleation. 

Frazil  forms  over  relatively  compact  length  and 
time  scales.  In  fact,  the  formation  phase  can  be  du¬ 
plicated  in  a  laboratory  beaker,  and  so  is  amenable 
to  laboratory  study.  Other  laboratory  experiments 
have  examined  the  rise  velocity  of  frazil  disks  and 
the  interaction  of  frazil  disks  with  suspended  sedi¬ 
ment.  Laboratory  work  has  also  been  inspired  by 
the  concept  of  desalination  of  sea  water  by  freezing. 
In  this  context,  frazil  ice  is  viewed  as  an  industrial 
product  that  is  manufactured  in  bulk  crystallizers. 

In  this  chapter,  we  will  concentrate  on  experi¬ 
ments  simulating  frazil  ice  in  natural  water  bodies. 

FRAZIL  ICE  FORMATION 

As  Chapter  4  describes,  the  basic  environmental 
parameters  controlling  frazil  ice  dynamics  are  the 
rates  of  heat  loss,  crystal  seeding  and  turbulent  en¬ 


ergy  dissipation.  Every  frazil  ice  experiment  has  re¬ 
quired  that  the  investigator  set  these  parameters 
and  then  observe  the  result.  Unfortunately,  investi¬ 
gators  often  have  not  been  able  to  directly  control 
the  seeding  rate,  but  rather  have  allowed  the  envi¬ 
ronmental  conditions  created  in  the  coldroom  to  do 
it.  So,  seldom  has  the  resultant  seeding  rate  been 
measured. 

The  heat  loss  rate  and  turbulent  energy  dissipa¬ 
tion  rate  can  be  determined  from  the  description  of 
the  experimental  apparatus  and  the  experimental 
procedure,  even  if  these  values  have  not  often  been 
reported.  These  last  two  parameters,  then,  provide 
us  with  a  straightforward  and  relatively  easy  way  of 
comparing  the  various  experiments  (see  Daly  and 
Axelson  [1989]  for  a  more  complete  analysis)  and  it 
is  usually  possible  to  calculate  these  values  as  de¬ 
scribed  below. 

The  heat  loss  rate  Q  may  be  determined  from  the 
rate  of  temperature  decline  in  the  turbulent  water 
before  frazil  forms.  To  determine  the  turbulent  en¬ 
ergy  dissipation  rate,  the  Reynolds  number  of  the 
flow  in  the  experiment  is  computed  from  the  aver¬ 
age  velocity  and  a  characteristic  length.  For  pipe 
flow,  the  diameter  of  the  pipe  is  the  characteristic 
length,  and  for  open-channel  flow,  the  hydraulic  ra¬ 
dius  is  the  characteristic  length.  The  friction  factor  is 
determined  from  the  Reynolds  number  of  the  flow 
using  an  estimated  surface  roughness.  For  pipe 
flow,  the  turbulent  energy  dissipation  rate  may  be 
found  using 


where  /  =  friction  factor 

u  =  average  velocity 
d  =  pipe  diameter. 

For  open-channel  flow,  the  turbulent  energy  dissi¬ 
pation  rate  is  a  function  of  the  friction  velocity  w», 
which  is  the  product  of  the  friction  factor  and  the 
mean  velocity 


(8) 
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The  turbulent  energy  dissipation  rate  is  then 


where  k  is  von  Karman's  constant,  generally  set 
equal  to  0.4,  and  Rh  is  the  hydraulic  radius. 

Table  1  contains  a  summary  of  data  from  the 
frazil  ice  experiments  described  below,  in  which  the 
water  temperature  record  as  a  function  of  time  was 
the  primary  result.  Frazil  ice  production  is  generally 
inferred  from  this  result,  a  reflection  of  the  great  dif¬ 
ficulty  in  measuring  the  concentration  or  size  distri¬ 
bution  of  the  frazil  ice  crystals  themselves.  Figure  2 
depicts  the  range  of  values  for  heat  loss  rate  and  tur¬ 
bulent  energy  dissipation  rate  found  from  the  ex¬ 
perimental  data. 

Michel  (1963)  carried  out  over  80  frazil  ice 
growth  experiments  in  an  outdoor  recirculating 
flume,  constructed  of  a  Plexiglas  channel,  30.5  cm 
wide  by  30.5  cm  deep  by  6.7  m  in  length,  connected 
to  a  15-cm-diameter  Plexiglas  pipe,  although  data 
were  reported  for  only  one  test.  A  20-cm-diameter, 
variable-speed  pump  with  a  cast  iron  casing  circu¬ 
lated  the  water  that  was  pumped  into  the  flume. 
Many  of  the  tests  were  done  at  night,  with  the  air 
temperature  ranging  from  -6.7  to  -32°C.  Water  tem¬ 
perature  was  measured  every  15  seconds  by  a  differ¬ 
ential  thermometer  that  was  accurate  to  ±0.0025°C. 


The  water  temperature  decline  rate  was  reported 
and  we  calculated  the  heat  loss  rate  from  this.  As¬ 
suming  that  the  experiment  was  dominated  by  the 
turbulence  immediately  downstream  of  the  pump, 
we  calculated  the  turbulent  energy  dissipation  rate 
for  this  flow  area.  Michel  reported  an  average  flow 
of  375  gal. /min  (0.025  m3/s),  yielding  a  velocity  in 
the  pipe  of  73  cm/s,  a  Reynolds  number  of  8.3X104 
and  a  friction  factor  of  0.023.  We  then  found  the  tur¬ 
bulent  energy  dissipation  rate  using  eq  7,  and  we 
measured  the  time  to  maximum  supercooling  from 
his  figure. 

Carstens  (1966)  studied  frazil  ice  in  a  recirculat¬ 
ing  oval  flume,  20  cm  wide  by  30  cm  deep  by  600  cm 
long,  located  in  a  -10°C  coldroom.  The  water  was  20 
cm  deep  and  was  circulated  by  a  variable  speed  pro¬ 
peller.  The  bottom  and  sides  of  the  acrylic  flume 
were  insulated,  and  cooling  of  the  water  was  aided 
by  a  fan  blowing  along  the  straight  portion  of  the 
flume.  Water  temperature  was  measured  with  a 
hand-held  mercury  thermometer  marked  to  0.01  °C 
that  was  immersed  from  5  to  10  cm.  In  a  typical  test, 
Carstens  measured  the  water  temperature  decline 
rate  and  then  determined  the  heat  loss  rate.  For  the 
experiments  reported,  the  average  flow  velocity  was 
50  cm/s,  yielding  a  Reynolds  number  of  1.8X104  and 
friction  factor  of  0.019.  We  calculated  the  turbulent 
energy  dissipation  rate  using  eq  9.  Unfortunately, 
for  several  of  the  experiments  shown  in  Carstens' 
Figure  5a,  the  turbulent  energy  dissipation  rate 
could  not  be  calculated  because  the  velocity  of  the 


Table  1.  Summary  of  experimental  data. 


t,  Q  e 


Experimenter 

Description  of  experiment 

Test  no. 

‘S 

(min) 

( j/cm 3  s) 

(cmVs1) 

Michel  (1963) 

Outdoor  recirculating  flume. 

27 

3.3 

1.2xl0-3 

33 

Carstens  (1966) 

Recirculation  flume  in  -10°C 

6A 

4 

1.9x1  tr3 

13 

coldroom;  propeller  in  tank  in 
-10°C  coldroom. 

5A 

8 

3.9x1  cr* 

Hanley  and 

Cylindrical  tank  with  paddles 

-2°C 

38.1 

0.9x10“* 

3 

Michel  (1977) 

in  coldroom.  Varied  coldroom 

-5°C 

21.6 

1.9X10-4 

3 

temperature  and  velocity. 

-10°C 

-20°C 

14.7 

8.1 

3x1CH 

7.1x10“* 

3 

Tsang  and  Hanley 

Warm  air  jacket  around  tank. 

A2(48%o) 

118.8 

1. 9xl0-3 

2 

(1985) 

except  test  C,  which  was  in 

Bl-3(23%») 

43.6 

2xl0-3 

2 

recirculating  flume.  Varied 

Bl-3(ll%o) 

36.4 

2xl0-3 

2 

salinity.  Seeded  with  shavings 

B3-4(fresh) 

18.1 

1.9x10-3 

2 

of  ice.  Average  coldroom 
temperature  -10°C. 

C-4(ocean) 

123.9 

1.1x10-3 

0.7 

Mueller  (1978) 

Supercooled  water  first  in 

E08 

— 

— 

1375 

agitation  tank  in  warm  water. 

E09 

■' '  f 

— 

4667 

6 


Hanley  &  Michel  (1977) 

Tsang  &  Hanley  (1985)  Tests  A-2 
B1-3 
B2-3 
B3-4 


Tsang  &  Hanley  (1985)  Test  C-4 


Carstens  (1 966)  Fig.  3-1 


Carstens  (1966)  Fig.  3-2 
Tsang  &  Hanley  (1985)  Tests  B1-3  &  B2-J 
Tsang  &  Hanley  (1985)  Tests  A7  &  B3-4 
Carstens  Fig.  3-3 

Michel  (1963),  Carstens(1966)  Test  5A 
Tsang  &  Hanley  (1985)  Test  C-4 
Carstens  (1966)  Fig.  3-4 
Hanley  &  Michel  (1977)  Air  T  -20”  C 
Carstens  (1966)  Fig.  3-5 

Carstens  (1966)  Test  6A 
Hanley  &  Michel  (1977)  Air  T  -10"  C 

Hanley  &  Michel  (1 977)  Air  T  -5*  C 


Hanley  &  Michel  (1977)  Air  T  -2”  C 


Figure  2.  Range  of  values  for  heat  loss  rate  and  turbulent  energy  dissipation  rate  from  summary  of  the 
frazil  ice  experiments'  data. 


water  was  not  reported.  However,  the  time  to  maxi¬ 
mum  supercooling  could  be  determined  from  this 
figure  and  is  included  for  comparison. 

Hanley  and  Michel  (1977)  conducted  frazil  ice 
experiments  in  a  stainless  steel  tank,  120  cm  in  diam¬ 
eter  and  76  cm  deep,  at  air  temperatures  of  -2,  -5, 
-10  and  -20°C.  They  did  not  report  the  water  depth 
in  the  tank.  Paddles  attached  to  an  axis  and  located 
25  cm  above  the  bottom  of  the  tank  circulated  the 


water  by  inducing  rotational  velocity.  For  these  ex¬ 
periments,  the  coldroom  temperature  was  set,  and 
then  the  paddles  started.  Water  temperature  was 
measured  by  a  thermistor  placed  just  below  the  wa¬ 
ter  surface.  We  calculated  the  heat  loss  rates  from  the 
reported  water  temperature  decline  rates.  The  re¬ 
sults  presented  by  Hanley  and  Michel  were  aver¬ 
aged  over  all  velocities;  the  minimum  velocity  at 
which  frazil  growth  was  noted  was  24  cm/s,  and  the 
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maximum  reported  velocity  was  73  cm/ s.  Assum¬ 
ing  a  water  depth  in  the  tank  of  70  cm  and  an  aver¬ 
age  velocity  of  48.5  cm/s,  we  calculated  the  Rey¬ 
nolds  number  at  8.8X104  and  estimated  the  friction 
factor  to  be  0.0145.  Equations  7  or  9  yielded  the  same 
turbulent  energy  dissipation  rate.  lime  to  maxi¬ 
mum  supercooling  was  not  reported;  total  time  of 
supercooling  was  reported,  however,  as  was  the 
time  from  the  minimum  temperature  to  the  end  of 
supercooling.  We  could  estimate  time  to  maximum 
supercooling  by  subtracting  the  latter  from  the 
former. 

Frazil  ice  formation  in  fresh  water,  ocean  water 
and  artificial  sea  water  of  varying  salinity  was  ex¬ 
amined  by  Tsang  and  Hanley  (1985).  They  placed  a 
rectangular  Plexiglas  tank,  38  cm  long  by  25.5  cm 
wide  by  15  cm  deep,  equipped  with  a  stirrer  at  one 
end  to  provide  turbulence,  in  a  -15°C  coldroom.  A 
horizontal  plate,  0.75  times  the  length  of  the  tank, 
was  placed  at  mid-depth  to  produce  vertical  recircu¬ 
lation  currents,  and  a  jacket  of  air,  slightly  above 
0°C,  was  placed  around  the  bottom  and  sides  of  the 
tank.  Cooling  was  provided  by  a  fan  that  pushed  air 
across  the  surface  of  the  water.  Tsang  and  Hanley 
measured  water  temperature  with  a  thermometer 
calibrated  to  0.0001°C,  with  repeatability  to  0.001°C. 
In  these  experiments,  the  water  was  first  cooled  to 
the  previously  selected  temperature  and  then  seed¬ 
ed  either  by  a  ball  of  ice  or  by  scrapings  from  ice. 

We  picked  one  sample  from  each  of  their  five 
groups  of  experiments  for  analysis  here.  We  calcu¬ 
lated  heat  loss  rates  for  each  sample  as  we  did  for 
Hanley  and  Michel.  All  but  the  group  C  (ocean  wa¬ 
ter)  experiments  took  place  in  the  tank  described 
above,  with  an  average  velocity  of  15  cm/ s.  Water 
depth  was  not  reported  and  we  assumed  it  to  be  15 
cm.  For  these  experiments,  we  calculated  the  Rey¬ 
nolds  number  to  be  3.9xl03,  with  a  resulting  friction 
factor  of  0.031,  and  we  estimated  turbulent  energy 
dissipation  rates  using  eq  7.  Group  C  tests  took 
place  in  a  recirculating  flume,  15  cm  wide  by  13  cm 
deep,  with  a  water  depth  of  11  cm,  and  average  ve¬ 
locity  of  15  cm/s.  Tsang  and  Hanley  reported  the 
Reynolds  number  as  8.54xl03.  From  this,  we  esti¬ 
mated  the  friction  factor  to  be  0.025,  and  we  esti¬ 
mated  the  turbulent  energy  dissipation  rate  using  eq 
9.  Time  to  maximum  supercooling  was  reported. 

Mueller  (1978)  studied  the  frazil  nucleation  pro¬ 
cess  using  an  agitating  tank,  17.2  cm  long  by  12  cm 
wide  by  20  cm  deep,  in  a  coldroom  kept  slightly 
above  0°C.  The  tank  was  surrounded  by  a  jacket 
through  which  a  coolant  was  circulated.  Agitation 
was  provided  by  a  grid  submerged  in  the  tank.  In 


these  experiments,  the  water  was  supercooled  to  the 
desired  degree,  the  coolant  circulation  was  stopped, 
agitation  was  begun  and  the  supercooled  water 
seeded.  Water  depth  ranged  from  20  to  20.5  cm.  He 
recorded  the  temperature  using  linear  type  ther¬ 
mistors  that  were  accurate  to  ±0.02°C  in  air  and 
±0.002°C  in  water.  Turbulent  energy  dissipation 
rates  were  reported.  These  experiments  were  not 
designed  to  allow  the  time  to  maximum  supercool¬ 
ing  to  be  determined,  as  the  water  was  not  seeded 
until  the  maximum  supercooling  was  achieved.  We 
report  the  turbulent  energy  dissipation  rates  only 
for  comparison  with  the  other  experiments. 

OTHER  ASPECTS  OF 
FRAZIL  ICE 

Daly  and  Colbeck  (1986)  conducted  a  series  of  ex¬ 
periments  designed  to  measure  the  size  distribution 
of  frazil  ice  crystals  at  four  equally  spaced  locations 
along  a  refrigerated  hydraulic  flume.  Several  differ¬ 
ent  slopes  and  bottom  roughnesses  were  tested  in 
an  approximately  -12°C  coldroom.  A  crystal  imag¬ 
ing  system  was  used  to  obtain  photographs  of  frazil 
crystals  within  a  control  volume  at  about  mid¬ 
depth.  Unfortunately,  physical  constraints  imposed 
by  the  experimental  conditions  (e.g.,  flow  depth) 
and  the  crystal  imaging  system  limited  measure¬ 
ments  to  the  longitudinal  number  distribution  of 
frazil  ice  crystals.  Vertical  number  distribution  could 
not  be  measured. 

Ice  crystal  diameters  determined  from  the  pho¬ 
tographs  ranged  from  35  pm  to  0.5  mm,  with 
mean  aspect  (diameter  to  thickness)  ratios  of  6.37, 
7.45,  8.53  and  9.61  at  the  four  measurement  sta¬ 
tions,  from  upstream  to  downstream.  The  number 
concentration  of  frazil  ice  crystals  ranged  from 
0.17  to  0.982/cm3.  Although  the  mean  frazil  crys¬ 
tal  size  and  number  concentration  generally  in¬ 
creased  in  the  downstream  direction,  they  did  not 
increase  consistently  for  each  test.  The  production 
of  anchor  ice,  which  acted  as  a  sink  for  suspended 
crystals,  was  thought  to  be  the  major  cause  for  the 
variation. 

Experiments  examining  the  distribution  or  con¬ 
centration  of  frazil  ice  in  the  water  column  have 
been  hindered  by  the  lack  of  accurate  measurement 
techniques.  In  fact,  in  a  review  of  frazil  ice  forma¬ 
tion,  Ostercamp  (1978)  stated  that  "..the  most  critical 
need  is  for  instrumentation  that  can  be  used  to  mea¬ 
sure  the  frazil  ice  concentration  in  a  cross  section  of  a 
river...".  A  number  of  researchers  (e.g.,  Kristinsson 
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1970)  have  used  changes  in  conductivity  resulting 
from  ice  formation  to  estimate  concentration;  how¬ 
ever,  their  results  appear  to  be  inconclusive.  A  calo¬ 
rimetric  method  of  measuring  frazil  ice  concentra¬ 
tion  currently  being  developed  has  been  tested  in 
the  laboratory  (Lever  et  al.  1992),  with  encouraging 
results. 

Ettema  et  al.  (1984)  examined  frazil  ice  nucle- 
ation  and  formation  in  a  rectangular  turbulence  jar 
equipped  with  an  oscillating  grid.  The  water  was 
supercooled  by  chiller  elements  containing  an  ethyl¬ 
ene  glycol  solution  at  -2°C,  and  natural  seeding  was 
limited  by  keeping  the  room  temperature  above 
freezing.  The  supercooled  water  was  artificially 
seeding  by  the  addition  of  a  small  block  of  ice  for 
most  experiments.  In  some  experiments,  a  syringe 
of  water  was  withdrawn  from  the  turbulence  jar,  ex¬ 
posed  to  a  piece  of  solid  ice,  and  then  returned  to  the 
jar,  presumably  containing  ice  nuclei. 

Despite  the  numerous  difficulties  typical  of  labo¬ 
ratory  frazil  ice  experiments,  the  researchers  were 
able  to  obtain  data  for  calibrating  their  analytical 
model.  Ettema  et  al.  reported  that  the  rate  of  frazil 
ice  formation  was  directly  proportional  to  the  turbu¬ 
lence  intensity  of  the  water  body.  Frazil  crystal  di¬ 
ameters  ranged  between  about  2  and  16  mm;  the 
larger  crystal  diameters  were  associated  with  in¬ 
creased  turbulence  intensity  and  also  with  greater 
initial  supercooling. 

The  rise  pattern  and  velocity  of  frazil  ice  crystals, 
parameters  important  in  frazil  entrainment,  trans¬ 
port  and  deposition,  have  been  observed  using 
frazil  ice  obtained  in  the  field  and  in  the  laboratory. 
Gosink  and  Ostercamp  (1983)  reported  on  a  series  of 
tests  conducted  in  the  Chatanika  River,  where  they 
collected  a  sample  of  river  water  in  a  transparent 
graduated  cylinder  and  then  observed  the  move¬ 
ment  of  the  visible  frazil  ice  crystals.  Rise  velocity, 
measured  by  timing  the  vertical  displacement  of  the 
crystals,  ranged  from  3  to  22  mm/s.  Crystal  diam¬ 
eters,  estimated  using  the  gradations  on  the  cylin¬ 
der,  ranged  from  1  to  6  mm. 

Wuebben  (1984)  used  a  similar  method  to  mea¬ 
sure  rise  velocities  of  artificial  frazil  ice  crystals  and 
ice  crystals  produced  in  a  hydraulic  flume  located  in 
a  coldroom.  He  measured  rise  velocity  against  the 
vertical  displacements  of  frazil  crystals  in  a  plastic 
cylinder,  finding  rise  velocities  of  about  0.8  to  5.4 
mm/s  for  estimated  0.4-  to  4-mm  crystal  diameters. 

The  interaction  between  rising  frazil  ice  crystals 
and  suspended  sediment  was  the  subject  of  a  study 
by  Kempema  et  al.  (1986).  They  did  both  freshwater 
and  saline  (0, 29  and  36  ppt  salinity)  tests;  this  chap¬ 


ter  discusses  the  freshwater  tests.  They  used  an  oval 
recirculating  flume,  located  in  a  laboratory  cold- 
room,  whose  sides  and  bottom  were  insulated  so 
that  heat  would  be  lost  primarily  through  the  sur¬ 
face  of  the  water.  A  4-cm  thick  layer  of  sand  was 
placed  on  the  bottom  of  the  flume,  where  the  water 
was  then  17  cm  deep.  Both  bed  load  and  suspended 
load  were  observed  under  the  range  of  induced  ve¬ 
locities  (42  to  70  cm/ s). 

Typical  frazil  disks,  1  to  5  mm  in  diameter,  joined 
together  to  form  roughly  spherical  floes  up  to  80 
mm  in  diameter.  These  floes  rose  toward  the  water 
surface,  but  some  touched  the  sand  bed  of  the  flume 
and  remained  there.  In  addition,  sediment  got 
trapped  within  the  floes,  at  times  in  quantities  suffi¬ 
cient  to  "sink"  them.  Kempema  et  al.  believed  that 
sediment  was  trapped  in  the  voids  between  the 
crystals,  rather  than  being  incorporated  into  the 
crystals  themselves. 

Suspended  sediment  samples  were  taken  before 
frazil  production  began  and  at  the  point  of  maxi¬ 
mum  frazil  production:  concentrations  ranged  from 
0.133  to  0.8  g/L  before  frazil  production,  and  from 
0.130  to  0.4  g/L  during  frazil  formation.  The  mean 
decrease  in  suspended  sediment  was  about  34% 
during  frazil  ice  production,  and  the  concentration 
of  sediment  contained  within  the  frazil  floes  ranged 
from  0.94  to  20.2  g/L. 

Kempema  et  al.  also  observed  the  interaction  be¬ 
tween  anchor  ice  and  sediment  during  these  experi¬ 
ments,  reporting  two  mechanisms  for  anchor  ice  for¬ 
mation:  frazil  floes  that  attached  to  objects  projecting 
into  the  flow  and  sediment-laden  frazil  floes  that 
settled  to  the  bottom.  Floes  that  settled  in  the  lee  of  a 
sand  ripple  were  soon  incorporated  into  the  bed 
structure  as  the  bedforms  migrated  in  an  ice-sedi¬ 
ment  block.  This  block  reappeared  if  the  ripples  mi¬ 
grated  past  the  buried  anchor  ice.  Sediment  concen¬ 
trations  measured  in  anchor  ice  samples  ranged 
from  4.74  to  37  g/L  and  from  17.6  to  42.8  g/L  in 
combination  frazil-anchor  ice  samples.  Again,  the 
authors  believed  that  sediment  was  trapped  within 
the  voids  between  the  crystals,  rather  than  incorpo¬ 
rated  into  the  crystals. 

Additional  experiments  using  the  same  appara¬ 
tus  were  reported  by  Kempema  et  al.  (1993)  for  both 
freshwater  and  saline  conditions.  As  before,  both 
suspended  and  bed  sediments  were  incorporated 
into  frazil  and  anchor  ice.  Concentrations  of  sus¬ 
pended  particulate  matter  ranged  from  0.94  to  20.2 
g/L  with  a  mean  of  8.59  g/L  in  frazil  floes  and  from 
4.74  to  42.8  g/L  with  a  mean  of  26.7  g/L  in  anchor 
ice  floes. 
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SUMMARY 

Laboratory  experiments  have  been  a  valuable 
tool  for  providing  insight  into  many  aspects  of  frazil 
ice  formation,  transport  and  interaction  with  sedi¬ 
ment.  The  lack  of  practical  and  efficient  ways  to 
measure  frazil  ice  concentrations  in  the  laboratory 
or  the  field  has  been  a  problem.  The  even  more  de¬ 


sirable  prospect  of  measuring  the  size  distribution 
of  frazil  crystals  seems  a  distant  goal.  The  small  size, 
nonspherical  shape,  and  optical  properties  of  frazil 
ice  present  daunting  obstacles  to  be  overcome  in 
pursuit  of  this  goal.  However,  direct  and  unambigu¬ 
ous  measurements  such  as  these  are  necessary  for 
progress  in  the  development  of  theoretical  models 
of  frazil  ice  to  continue. 
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CHAPTER  3 

Evolution  of  Frazil  Ice  in  Natural  Water  Bodies 


STEVEN  F.  DALY 


INTRODUCTION 

One  chief  characteristic  of  the  "peculiar  mode  of 
ice  formation"  in  natural  water  bodies  known  as 
frazil  ice  is  the  continuous  evolution  of  its  form.  The 
major  paths  of  this  evolution  are  diagrammed  in 
Figure  3.  Arguably,  there  are  many  different  ways  of 
describing  this  evolution,  but  for  the  sake  of  clarity, 
brevity  and  completeness,  we  will  identify  three 
general  phases  of  frazil  evolution. 

The  first  phase  is  formation,  characterized  by  su¬ 
percooled  water,  turbulent  flow,  the  rapid  growth  of 
disk-shaped  crystals,  and  the  creation  of  new  crys¬ 
tals  by  secondary  nucleation.  The  length  scales  of 
the  ice  associated  with  this  stage  (Fig.  4)  range  from 
several  micrometers  to  perhaps  a  few  millimeters. 
This  stage  usually  lasts  for  a  relatively  short  time 
during  very  cold  periods  when  the  heat  transfer 
from  open  water  surface  is  large. 

The  second  phase  is  transformation  and  trans¬ 
port.  This  phase  follows  the  first  in  time,  and  results 
from  the  rapid  production  of  the  first  phase.  It  is 
characterized  by  water  more  or  less  at  the  equilib¬ 
rium  temperature,  and  frazil  in  the  form  of  floes,  an¬ 
chor  ice  and  floes.  The  length  scales  of  the  ice  associ¬ 
ated  with  this  stage  (Fig.  4)  range  from  several  milli¬ 
meters  to  several  meters.  The  frazil  is  largely  mov¬ 


ing  under  the  influence  of  the  river  or  stream,  gener¬ 
ally  at  the  surface.  This  ice  may  travel  long  distances 
and  move  for  many  days.  After  cold  nights,  it  is  typi¬ 
cal  to  see  slush,  formed  of  frazil  floes,  moving  along 
at  the  water  surface  of  northern  rivers  and  streams. 
This  slush  may  eventually  form  large  moving  floes. 

The  third  phase  is  characterized  by  stationary, 
floating  ice  covers  that  may  be  quite  large  and  last 
for  the  entire  winter  season.  These  ice  covers  are 
formed  by  a  variety  of  mechanisms,  depending  on 
the  form  of  the  frazil  ice  when  it  arrives  at  the  sta¬ 
tionary  ice  cover  and  the  hydraulic  conditions  at  the 
cover's  leading  edge.  These  floating  covers  may 
raise  stages  and  cause  flooding,  cause  increases  in 
the  channel  head  losses  that  can  disrupt  power  pro¬ 
duction,  and  interfere  with  navigation. 

To  focus  our  discussion,  we  can  further  subdi¬ 
vide  the  three  phases  discussed  above  into  six 
process  categories  (Fig.  3).  Each  category  is  distin¬ 
guished  by  one  or  two  distinct  forms  of  frazil  ice 
and  one  or  two  dominating  processes.  The  pro¬ 
gression  of  the  categories  in  general  follows  the 
known  evolution  of  frazil  ice.  We  can  outline  each 
category  as  follows. 

In  the  formation  phase  there  are  two  categories. 
Seeding  describes  the  introduction  of  ice  crystals 
from  outside  the  water  body  that  starts  the  forma- 
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Figure  3.  Evolution  of  frazil  ice  in  natural  water  bodies. 
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tion  of  frazil  ice.  Frazil  ice  dynamics  describes  the 
interaction  of  seeding,  secondary  nucleation  and 
crystal  growth.  It  is  during  this  stage  that  disks  ap¬ 
pear,  each  being  a  separate  crystal.  Disk  crystals  are 
the  unique  form  of  frazil  ice  and  the  form  by  which 
it  is  chiefly  known.  These  disks  are  formed  in  super¬ 
cooled,  turbulent  water.  A  complete  discussion  of  the 
physics  and  dynamics  of  this  stage  can  be  found  in 
Chapters  1  and  4. 

In  the  transformation  and  transport  phase,  frazil 
ice  displays  a  marked  tendency  to  aggregate  into 
larger  and  larger  masses.  At  this  time,  we  cannot  de¬ 
termine  the  aggregation  rate,  or  even  if  the  progres¬ 
sion  to  larger  sizes  will  happen  at  all.  In  general,  the 
progression  goes  at  a  rate  determined  by  the  turbu¬ 
lent  intensity  of  the  flow,  but  very  high  levels  of  tur¬ 
bulence  will  keep  the  larger  forms  of  frazil  from  oc¬ 
curring.  This  means  that  the  particular  evolutionary 
pathway  followed  by  the  frazil  in  its  progression  to 
larger  aggradations  will  be  unique  for  each  water 
body.  As  a  result,  it  is  quite  difficult  to  generalize  the 
transformation  and  transport  stage  of  frazil  evolu¬ 
tion.  The  three  categories  that  occur  in  this  stage — 
flocculation  and  deposition,  transport  and  mixing, 
and  floe  formation  and  induration — describe  the 
processes  through  which  the  initial  disk  crystals  are 
transformed  into  larger  and  larger  ice  masses.  The 
length  scales  associated  with  these  forms  can  be¬ 
come  quite  large  (Fig.  4). 

The  small  disk  crystals  of  frazil  ice,  which  start  at 
the  size  of  millimeters,  can  eventually  form  very 


large  river  ice  covers  with  lengths  of  many  kilome¬ 
ters  and  thicknesses  of  many  meters  (the  stationary 
ice  cover,  phase  six). 

We  will  now  discuss  each  of  the  categories  and 
the  various  forms  of  frazil  that  define  each. 

SEEDING 

For  a  long  time,  the  origin  of  frazil  ice  under 
natural  conditions  was  debated.  Spontaneous  nu¬ 
cleation  of  ice,  either  heterogeneously  or  homoge¬ 
neously,  was  thought  responsible  for  the  ultimate 
origin  of  frazil.  All  available  data  indicate  that  spon¬ 
taneous  nucleation  of  ice  is  not  possible  in  natural 
water  bodies;  therefore,  seed  crystals  are  necessary 
(Daly  1984).  Unfortunately,  there  is  very  little  quan¬ 
titative  information  available  on  either  the  sizes  of 
the  seed  crystals  or  the  seeding  rates.  Osterkamp 
(1977)  reports  that  hexagonal  plate  ice  crystals,  rang¬ 
ing  from  60  to  350  pm,  were  observed  in  the  air 
above  a  supercooled  Alaskan  stream.  Their  concen¬ 
trations  ranged  from  6  to  6X104  crystals/m3.  It  is 
possible  to  estimate  a  surface  seeding  rate  from  the 
concentration  of  crystals  suspended  in  the  air  CA  as 

4  =  CA«f  (1) 

where  u(  is  the  terminal  fall  velocity  of  the  ice  crys¬ 
tals.  The  terminal  fall  velocities  of  the  largest  and 
smallest  ice  crystal  are  estimated  as  17  and  5  cm/ s, 
respectively  (Pruppacher  and  Klett  1978),  yielding  a 
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seeding  rate  per  unit  surface  area  of  3xl0“s  to  1.0 
crystals/cm2  s.  For  a  discussion  of  other  mechan¬ 
isms  for  producing  seed  crystals  see  Chapter  1. 

FRAZIL  ICE  DYNAMICS- 
DISK  CRYSTALS 

We  know  from  observations  that  the  dominant 
shape  of  ice  crystals  that  grow  at  the  supercooling 
levels  found  in  turbulent  water  bodies  is  a  flat  disk. 
Virtually  all  field  observations  of  frazil  ice  note  that 
the  crystals  are  disk  shaped.  It  has  been  reported 
that  ice  crystals  in  the  shape  of  six-pointed  stars, 
hexagonal  plates  or  spheres,  and  small  pieces  of 
dendritic  ice  all  evolve  into  the  disk  shape  in  natural 
water  bodies. 

Disk-shaped  crystals  have  been  studied  in  the 
laboratory  by  Kumai  and  Itagaki  (1953),  Arakawa 
(1954)  and  Williamson  and  Chalmers  (1966).  Ara¬ 
kawa  created  disk  crystals  by  first  growing  den¬ 
dritic  ice  crystals  on  the  bottom  of  a  container 
containing  slightly  supercooled  water.  He  then 
scratched  the  crystals  with  the  tip  of  a  fine  needle. 
Spherical  ice  particles  with  a  diameter  of  about  IQ-2 
mm  rose  towards  the  surface;  two  flat  spots  formed 
on  the  spherical  particles'  surface  as  they  floated  up. 
Further  growth  was  lateral;  the  particles  became 
disks.  The  crystals  grew  into  disks  with  diameters  of 
0.5  to  3  mm  and  with  diameter-to-thickness  ratios  of 
5:1  to  100:1.  Once  the  disks  grew  to  a  certain  size 
they  took  on  a  notched  look,  with  many  notch¬ 
shaped  growths  at  their  perimeters. 

The  disk  shape  of  the  ice  crystal  is  morphologi¬ 
cally  instable.  In  natural  water  bodies,  this  insta¬ 
bility  is  observed  as  scalloped  edges,  dendritic 
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growths,  irregular  protuberances,  etc.,  on  the  perim¬ 
eter  of  the  crystals.  Experiments  have  shown  that 
the  disk  shape  always  becomes  unstable  if  the  su¬ 
percooling  is  large  or  if  the  diameter  of  the  crystal  is 
large.  Instability  is  perhaps  the  mechanism  that  lim¬ 
its  the  maximum  disk  radius  size.  The  entire  reason 
for  the  instability  of  the  disk  shape  is  not  known, 
although  it  has  been  the  subject  of  several  investiga¬ 
tions  (Williamson  and  Chalmers  1966,  Fujioka  and 
Sekerka  1974).  Williamson  and  Chalmers  concluded 
that  the  instability  of  the  disk  shape  depends  on  heat 
flow  into  the  liquid,  not  crystallography.  This  con¬ 
clusion  seems  correct,  as  the  diameter  at  which  the 
crystal  becomes  unstable  varies  between  experi¬ 
ments. 

Disk-shaped  frazil  ice  crystals  have  been  ob¬ 
served  under  three  contexts:  in  natural  water  bod¬ 
ies,  in  laboratory  experiments  and  in  research  into 
desalination  by  freezing.  The  difficulties  in  measur¬ 
ing  frazil  size  or  concentration  have  limited  the 
number  of  measurements  made. 

Estimates  of  the  range  of  frazil  crystal  concentra¬ 
tions  in  natural  rivers  is  from  about  10,000  to  1  mil¬ 
lion  crystals/m3  (Schaefer  1950,  Osterkamp  and 
Gosink  1982).  Measurements  in  a  laboratory  flume 
(Daly  and  Colbeck  1986)  indicated  concentrations 
ranging  from  1.80x10s  to  9.82x10s  crystals/m3. 

The  sizes  of  crystals  have  been  measured  in  the 
range  of  0.05  mm  to  several  millimeters  (Daly  and 
Colbeck  1986).  The  aspect  ratio  of  the  crystals  is 
shown  in  Figure  5.  The  observations  of  Margolis 
(1969)  on  frazil  ice  grown  during  desalination  by 
freezing  indicate  that  the  thickness  of  the  frazil  disk 
is  (0.68  ±  0.17)R,  where  R  is  the  major  radius  of  the 
disk.  Smith  and  Sarofin  (1979)  report  that  the  maxi¬ 
mum  diameter  is  approximately  1.6  mm  for  disk 

crystals  produced  in  turbulent  crystallizers. 


Figure  5.  Size  distribution  of  the  major  crystal  diameters 
measured  in  a  laboratory  flume  (from  Daly  and  Colbeck  1986). 


FLOCCULATION 
AND  DEPOSITION 

Frazil  floes 

The  flocculation  of  frazil  ice  is  not  well  under¬ 
stood,  although  sintering  is  one  mechanism  used 
to  explain  it.  Sintering  apparently  results  from 
the  tendency  of  crystals  to  minimize  their  surface 
energy.  Martin  (1981)  stated  that  for  disks  of 
thickness  on  the  order  of  1-10  pm,  the  time  for 
bonds  to  form  is  fast,  perhaps  10~2  seconds. 
Mercier  (1984),  in  his  simulations,  assumed  that 
turbulent  shear  and  differential  rising  were  the 
two  mechanisms  causing  collisions  of  crystals  for 
sintering. 
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Descriptions  of  floes  in  streams  (Osterkamp  and 
Gosink  1982)  and  observations  by  this  author  sug¬ 
gest  that  floes  are  formed  initially  by  a  few  to  several 
hundred  frazil  disks.  The  bonds  within  newly 
formed  floes  are  initially  very  weak,  and  they  will 
often  collapse  when  removed  from  the  water.  Large 
floes,  called  frazil  pebbles,  were  described  by 
Chacho  et  al.  (1986)  as  ranging  in  length  from  less 
than  5  to  greater  than  150  mm.  Grains  composing 
the  frazil  pebbles  do  not  possess  a  preferred  c-axis 
orientation,  but  appear  to  show  alignment  related  to 
grain  size  and  shape.  This  may  imply  that  these 
pebbles  are  not  formed  by  collisions  caused  by  tur¬ 
bulent  shear,  which  would  presumably  result  in  ran¬ 
domly  oriented  grains,  but  probably  result  from  dif¬ 
ferential  rising  or  a  depositional  mechanism,  such  as 
might  take  place  during  the  formation  of  anchor  ice. 
The  round  shape  of  the  pebbles  may  result  from 
mechanical  abrasion,  indicative  of  the  long  distance 
that  they  may  have  traveled  on  the  Tanana  River  be¬ 
fore  they  were  observed,  or  represent  the  equilib¬ 
rium  form  of  floes  that  have  traveled  long  distances 
under  water. 

Anchor  ice 

There  has  been  little  quantitative  study  (Mar- 
cotte  and  Robert  1986)  of  anchor  ice,  which  is 
formed  initially  by  the  deposition  of  suspended 
frazil  crystals.  Once  deposited,  the  heat  transfer  rate 
from  the  anchor  ice  is  much  greater  than  that  for  the 
suspended  frazil.  The  importance  of  anchor  ice  as  a 
sink  of  suspended  crystals  and  as  a  source  of  latent 
heat  in  an  overall  stream  crystal  balance  or  heat  bal¬ 
ance  will  be  inversely  proportional  to  the  hydraulic 
radius  of  the  stream.  Anchor  ice  may  be  the  domi¬ 
nant  form  of  frazil  produced  in  shallow  streams. 

TRANSPORT  AND  MIXING 
Frazil  slush 

Frazil  slush  is  the  collection  of  frazil  floes  and  in¬ 
dividual  frazil  crystals  on  the  water  surface  in  a  dis¬ 
tinct  layer.  The  mechanical  properties  of  the  highly 
porous  frazil  slush  are  not  known.  A  sample  of  frazil 
slush  lifted  from  the  water  will  quickly  drain,  re¬ 
vealing  a  characteristic  white  appearance.  On  the 
surface  it  has  a  marked  tendency  to  form  a  "struc¬ 
ture"  of  independent,  individual  clumps,  that  even¬ 
tually  become  floes. 

How  this  structure  arises  is,  at  present,  un¬ 
known.  It  may  well  arise  "spontaneously"  as  nu¬ 
merical  simulations  of  cohensionless,  inelastic  par¬ 
ticles  suggest  (Hopkins  and  Lange  1991),  or  it  may 


result  from  the  buffeting  of  the  slush  by  surface  ed¬ 
dies.  The  degree  to  which  the  formation  of  a  struc¬ 
ture  is  aided  by  the  mechanical  interlocking  of  the 
rough  frazil  floes  is  also  unknown.  In  any  case,  the 
initial  clumps  that  are  formed  can  be  broken  up 
upon  reaching  a  highly  turbulent  portion  of  river. 

Tsang  (1988)  reported  15%  concentrations  of  ice 
in  the  surface  layer  at  the  Cachine  Rapids  of  the  St. 
Lawrence  River.  Matousek  (1981)  cited  densities  of 
frazil  slush  from  300  to  640  kg/m3,  density  and  po¬ 
rosity  being  linearly  dependent  on  the  mean  air 
temperature. 

Suspension 

The  form  of  the  frazil  ice  that  is  in  suspension  can 
vary  so  widely  that  it  is  difficult  to  adequately  de¬ 
scribe  it.  The  form  will  chiefly  depend  on  the  frazil's 
past  history.  Close  to  the  point  of  formation,  more 
disk  crystals  will  be  found  in  the  suspension;  far 
from  the  point  of  formation,  there  will  be  more  floes. 
Another  consideration  is  the  relative  time  that  the 
ice  has  spent  on  the  surface  as  compared  to  being  in 
the  depth  of  the  flow.  If  the  frazil  has  traveled  on  the 
surface  and  then  been  resuspended  (by  passing  into 
a  zone  of  high  turbulence,  for  example),  a  very  high 
proportion  of  floes  could  be  expected. 

Arden  and  Wigle  (1972)  observed  frazil  ice  in 
suspension  on  the  Niagara  River,  where  frazil  ice 
was  largely  in  the  form  of  disks.  They  viewed  the 
suspension  with  an  underwater  lamp  at  night,  giv¬ 
ing  "the  viewer  the  impression  that  he  is  looking  at  a 
driving  snowstorm,  as  seen  through  the  headlight 
beams  of  a  moving  automobile  at  night." 

There  are  very  few  measurements  of  the  concen¬ 
tration  of  frazil  ice  in  suspension.  One  profile,  mea¬ 
sured  by  Tsang  (1986)  using  a  device  that  inferred 
the  concentration  of  frazil  based  on  the  electrical 
conductivity  of  the  water,  gave  values  in  the  range 
of  0  to  0.03%. 

Determination  of 
vertical  transport 

To  adequately  discuss  the  vertical  mixing  of 
frazil  ice  requires  that  the  size  distribution  and 
form  of  the  frazil  crystals  be  well  known,  that  the 
buoyancy  of  the  frazil  can  be  determined,  that  the 
turbulence  of  the  water  flow  field  be  well  de¬ 
scribed,  that  the  heat  transfer  rate  from  the  water 
surface  be  known,  and  that  the  evolution  in  form 
of  the  frazil  in  response  to  these  influences  can  be 
estimated.  Unfortunately,  at  present,  we  know 
very  little  of  the  above.  Because  of  this,  investiga¬ 
tors  have  concentrated  on  specific  questions,  such 
as:  What  is  the  buoyant  rise  velocity  of  frazil  disks 
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and  frazil  floes?  Under  what  conditions  of  flow 
and  heat  loss  can  a  distinct  surface  layer  of  slush 
be  expected  to  form?  What  is  the  concentration  of 
frazil  ice  in  suspension? 

Matousek  (1984)  determined  the  flow  condi¬ 
tions  under  which  a  distinct  surface  layer  could  be 
expected  by  comparing  the  buoyant  rise  velocity  of 
frazil  disks  and  the  mean  magnitude  of  the  turbu¬ 
lent  velocity  fluctuations  in  the  vertical  direction.  If 
the  rise  velocity  exceeded  this  mean  magnitude, 
then  a  surface  layer  should  form.  Empirical  correla¬ 
tions  were  used  to  determine  both  the  rise  velocity 
and  the  vertical  velocity  fluctuations. 

Gosink  and  Osterkamp  (1983)  compared  the 
buoyant  and  mixing  time  scales  to  determine  "well 
mixed  vs.  layered"  flow.  If  the  buoyant  time  scale 
(defined  by  the  time  required  for  a  disk  to  rise 
through  the  entire  depth)  was  less  than  the  time  re¬ 
quired  for  the  frazil  disk  to  diffuse  downward 
through  the  entire  depth,  then  they  expected  a  dis¬ 
tinct  surface  layer  to  form.  Interestingly,  the  ratio  of 
the  time  scales  is  independent  of  the  flow  depth,  and 
depends  only  on  the  mean  flow  velocity,  the  channel 
roughness  and  the  rise  velocity  of  the  frazil.  They 
assumed  the  rise  velocity  to  be  0.01  m/s,  which  is 
appropriate,  based  on  their  analysis,  for  disks  2.5 
mm  in  diameter,  with  a  thickness-to-diameter  ratio 
of  1:10. 

Tsang  (1988)  compared  the  kinetic  energy  avail¬ 
able  for  mixing  and  the  increase  in  potential  energy 
resulting  from  suspension  of  the  frazil  ice  in  the 
flow.  From  this  he  derived  a  "distribution  equa¬ 
tion,"  which,  when  solved  numerically,  provided 
the  vertical  concentration  distribution  of  the  frazil 
ice,  including  a  distinct  surface  layer,  if  present. 

Lai  (1989)  and  Gunaratna  (1989)  developed  "two 
layer"  models,  with  one  layer  corresponding  to  the 
distinct  surface  layer  and  the  other  to  the  suspended 
frazil.  Essentially,  the  ice  exchange  between  the  lay¬ 
ers  was  determined  by  one  parameter  whose  value 
was  found  through  model  calibration. 

Ho  (1990)  determined  the  vertical  diffusion  of  the 
frazil  through  the  water  depth  by  solving  the  Fick- 
ian  diffusion  equation.  Several  empirically  deter¬ 
mined  coefficients  were  required  in  this  approach. 

Liou  and  Ferrick  (1992)  developed  a  model  in 
which  the  net  upward  migration  of  frazil  disks  at¬ 
tributable  to  buoyancy  was  opposed  by  intermittent 
mixing  induced  by  large  energy-containing  eddies. 
A  surface  renewal  model  was  used  to  describe  the 
effects  of  large  eddy  mixing.  They  supposed  that  a 
critical  surface  layer  thickness  was  necessary  for  the 
surface  to  develop  structure.  In  general,  they  found 
that  the  heat  loss  rate  from  the  water  surface,  the 


surface  renewal  frequency,  and  the  critical  surface 
layer  thickness  determined  whether  the  frazil 
evolved  towards  a  well-mixed  equilibrium  state  or  a 
layered  state. 

FLOE  FORMATION 

How  floes  form  from  frazil  ice  floes  is  probably 
one  of  the  least  well  described  processes  in  the  evo¬ 
lution  of  frazil  ice.  As  mentioned  previously,  frazil 
slush  at  the  water  surface  shows  a  marked  tendency 
to  clump  together.  The  initial  clumps,  if  they  remain 
on  the  surface  long  enough,  freeze  together  and 
form  pans,  or  small  floes.  These  pans  often  grind 
against  one  another,  become  roughly  circular  and 
gain  upturned  edges,  and  then  are  known  as  pan¬ 
cake  ice.  The  pans  can  also  stick  together  and  form 
much  larger  floes,  which  attain  sizes  comparable  to 
the  river's  width. 

Pans  are  always  much  smaller  than  the  stream 
width,  and  as  a  result,  their  form  is  essentially  inde¬ 
pendent  of  the  stream  geometry.  Larger  floes,  on  the 
other  hand,  result  from  changes  in  the  stream  sur¬ 
face  velocities  that  tend  to  concentrate  pans  at 
specific  points  along  the  stream  channel.  The  form 
of  these  large  floes  then  is  strongly  controlled  by  the 
particular  flow  conditions  of  the  streams. 

The  formation  of  floes  of  all  sizes  is  governed  by 
two  separate  but  intertwined  processes:  the  increase 
in  floe  area  and  floe  induration.  The  increase  in  floe 
area  reflects  the  accumulation  of  frazil  ice  into  a 
single,  moving  unit.  Floe  induration  describes  the 
increase  in  floe  strength  because  of  exposure  to  the 
cold  air.  Immediately  after  their  formation,  floes  are 
quite  porous  and  have  little  strength.  Eventually, 
they  may  become  more  or  less  solid  ice  as  the  inter¬ 
stitial  water  freezes.  The  relation  among  the  floe  po¬ 
rosity,  size  and  strength  is  not  known  at  the  present. 
The  only  measurements  of  floe  strength  are  those  by 
White  (in  prep.),  who  estimated  in-situ  floe  strength 
by  suspending  instrumented  rods  in  a  mountain 
stream  during  a  period  when  frazil  ice  floes  were 
present.  She  observed  three  types  of  interactions  be¬ 
tween  the  rods  and  the  ice  floes:  a  "crushing-type" 
interaction,  in  which  the  rod  left  an  indentation 
track  in  the  floe;  a  liquid-like  interaction,  where  the 
floes  "flowed  around  the  rod"  and  resumed  their 
original  shape  with  no  discernible  tracks;  and  solid- 
impact  interactions,  where  the  floes  bounced  off  the 
rod.  She  was  able  to  draw  only  qualitative  infer¬ 
ences  between  the  type  of  interaction  and  the  air 
temperature,  travel  distance  and  appearance  of  the 
floes. 
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Osterkamp  and  Gosink  (1982)  described  floe  for¬ 
mation  in  interior  Alaskan  streams  in  detail.  They 
identified  five  mechanisms  that  can  produce  floes 
from  smaller  frazil  pans:  1)  contact,  penetration  and 
bonding  of  the  frazil  pans;  2)  compaction  and  drag 
cutoff*;  3)  compaction  and  convergent  flow  with 
cutoff  by  impact  of  incoming  pans;  4)  extrusive  flow 
and  drag  cutoff;  and  5)  agglomeration  with  cutoff 
controlled  primarily  by  river  curvature.  Extrusive 
flow  and  drag  cutoff  produced  very  large  floes  on 
the  Yukon  River,  with  maximum  dimensions  ex¬ 
ceeding  1  km. 

ICE  COVER  FORMATION 
AND  UNDER-ICE  TRANSPORT 

The  formation  of  river  ice  covers  by  frazil  ice  pro¬ 
ceeds  by  a  number  of  processes  (Ashton  1986, 
Michel  1971,  Osterkamp  and  Gosink  1982).  The  par¬ 
ticular  process  at  any  given  point  is  determined 
largely  by  the  form  of  the  frazil,  the  hydraulic  condi¬ 
tions  of  tire  flow,  and  the  heat  transfer  rate  from  the 
water  surface.  The  cover  can  be  considered  initiated 
when  the  frazil  ice  first  stops  because  of  an  ice  or 
hydraulic  control  structure,  because  of  the  growth  of 
surface  ice  in  a  slow  moving  reach,  or  because  the 
moving  frazil  ice  has  "arched"  across  the  open  wa¬ 
ter.  The  arching  of  floes  has  been  the  subject  of  some 
investigation,  but  arching  of  frazil  slush  has  re¬ 
ceived  much  less  attention,  although  its  occurrence 
in  the  field  has  been  reported  (Osterkamp  and 
Gosink  1982).  That  frazil  slush  can  arch  implies  that 
it  has  the  ability  to  transmit  shear  stresses,  but  the 
effective  viscosity  of  slush  and  its  other  material 
properties  are  not  known. 

A  major  issue  for  hydraulic  engineers  studying 
the  formation  of  ice  covers  is  determining  the  condi¬ 
tions  for  stability  of  the  frazil  ice  arriving  at  the 
cover's  leading  upstream  edge.  The  arriving  ice  is 
said  to  be  stable  if  it  retains  orientation  and  shape  at 
the  time  that  it  stops  at  the  leading  edge. 

The  form  of  frazil  ice  arriving  at  the  upstream 
end  of  an  ice  cover  plays  a  major  role  in  determining 
if  and  how  the  frazil  will  be  incorporated  into  it. 
Fra7.il  ice  that  is  largely  in  suspension  will  pass  un¬ 
der  the  cover  and  deposit  on  the  underside  at  some 
point  downstream.  Frazil  slush  may  be  incorpo- 


*  Cutoff  is  the  process  whereby  distinct  floes  are  sepa¬ 
rated  from  a  large  mass  of  compacted  frazil  ice.  Drag 
cutoff  refers  to  separation  that  occurs  when  the  fluid 
drag  acting  on  the  ice  exceeds  the  tensile  strength  of  the 
compacted  ice  mass. 


rated  into  the  ice  cover  at  the  surface  if  the  flow  ve¬ 
locity  is  low,  may  "pack"  at  higher  velocities,  or  may 
be  carried  under  the  cover  at  still  higher  velocities. 
Floes  may  be  "stable"  (that  is,  they  do  not  under- 
turn)  at  the  leading  edge  of  the  ice  cover  if  the  veloci¬ 
ties  are  low,  undertum  and  remain  at  the  leading 
edge  at  higher  velocities,  and  be  carried  under  the 
ice  cover  at  still  higher  velocities.  Determining  the 
range  of  hydraulic  conditions  at  which  each  of  these 
behaviors  will  occur  has  been  one  of  the  major  direc¬ 
tions  of  river  ice  research  (for  overviews,  see  for  ex¬ 
ample,  Michel  [1978]  or  Ashton  [1986]). 

Frazil  ice  that  is  transported  under  an  ice  cover 
often  deposits  on  the  underside,  and  there  is  a  long 
history  of  researchers  considering  this  frazil  ice 
deposition  to  be  analogous  to  that  of  sediment 
(Pariset  and  Hausser  1961,  etc.).  How  closely  this 
analogy  should  be  drawn  is  not  clear;  for  example, 
many  of  the  bed  forms  found  in  rivers  transporting 
sediment,  such  as  dunes,  antidunes,  ripples,  etc., 
have  not  been  reported  on  the  underside  of  river  ice 
covers  transporting  frazil.  (It  must  be  noted,  how¬ 
ever,  that  ripples  have  been  observed  to  form  in 
solid  ice  covers.  This  results  from  the  turbulent  flow 
beneath  the  cover  and  the  associated  heat  transfer 
[Ashton  1986],  not  frazil  transport.)  The  one  form 
that  has  been  observed,  which  can  grow  to  spectacu¬ 
lar  depths,  is  called  a  "hanging  dam"  (Michel  et  al. 
1981,  Beltaos  and  Dean  1981,  Hopper  et  al.  1981, 
Shen  et  al.  1984,  Michel  1990),  but  this  form  has  no 
recognized  analogy  in  sediment  transport. 

In  any  case,  applications  of  concepts  from  sedi¬ 
ment  transport  have  not  proven  especially  useful  in 
predicting  the  under-ice  deposition  of  frazil — other 
than  the  relatively  simple  and  straightforward  con¬ 
cept  of  a  "critical  velocity  criterion"  to  determine  the 
limiting  thickness  of  frazil  deposition  (Tesaker  1975, 
Michel  and  Drouin  1975,  Shen  et  al.  1984,  Michel 
1984,  Sun  and  Shen  1988).  Unfortunately,  the  value 
of  this  critical  velocity  has  been  found  to  vary  from 
about  0.9  m/s  at  the  beginning  of  the  winter  to 
about  0.5  m/s  at  the  end  (Michel  et  al.  1981)  at  one 
location  in  Canada,  and  to  vary  from  0.4  to  1.2  m/ s 
from  reach  to  reach  in  the  Yellow  River  in  China 
(Sun  and  Shen  1988).  Undoubtedly,  the  value  of  this 
critical  velocity  is  influenced  by  the  form  of  the 
frazil,  the  distance  it  has  traveled  under  the  ice 
cover,  the  degree  of  supercooling  of  the  river  and 
perhaps  other  factors. 

The  response  of  frazil  slush  on  arriving  at  the 
leading  edge  of  an  ice  cover  has  been  described  in 
several  works  of  B.  Michel  (1991,  1986,  1984),  in 
which  frontal  progression  and  packing  are  distin¬ 
guished  and  analyzed  through  the  use  of  momen- 
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turn,  energy  and  mass  flux  balances  at  the  leading 
edge.  The  stability  of  floes  arriving  at  the  leading 
edge  has  most  recently  been  analyzed  by  Daly  and 
Axelson  (1990)  and  McGilvary  and  Coutermarsh 
(1992). 

SUMMARY 

The  evolution  of  frazil  in  rivers  and  stream  is  in¬ 
fluenced  by  the  structure  and  growth  kinetics  of 


frazil  ice  crystals,  the  turbulent  intensity  of  the  flow, 
the  surface  heat  loss  rate  to  the  atmosphere,  the  plan 
form  of  the  river,  the  existence  of  hydraulic  and  ice 
control  structures,  and  many  other  factors.  In  broad 
outline,  the  evolutionary  pathway  of  frazil  is  well 
known,  but  little  of  it  can  be  quantitatively  de¬ 
scribed.  At  present,  efforts  to  control  or  modify  this 
evolution  must  rely  on  empirical  criteria  that  have 
been  obtained  through  experience  and  long  obser¬ 
vation. 
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CHAPTER  4 
Frazil  Ice  Dynamics 

STEVEN  F.  DALY 


NOMENCLATURE 

B  birth  function 

C  heat  capacity  of  fluid 

Q  heat  capacity  of  ice 

Cj  fluid  impurity  concentration 

D  death  function 

Erb  energy  created  by  crystal-boundary 

collisions 

E{v\,r2)  energy  created  by  collisions  of  crystals  of 
size  i\  and  r2 

Et  rate  of  energy  transfer  by  collision 

F|  number  of  particles  generated  per  unit  of 

collision  energy 

E2  fraction  of  particles  surviving  to  become 

crystals 
g '  gravity 

g'  reduced  gravity 

G  crystal  growth  rate 

h  heat  transfer  coefficient 

k  heat  conductivity  of  fluid 

Kv  crystal  shape  factor 

L  mean  latent  heat  of  fusion  of  ice 

m*  nondimensional  crystal  size 

m(r)  mass  of  crystal  of  size  r 

n  size  distribution  function 

N  total  number  of  crystals  per  unit  volume 

N  i  rate  of  introduction  of  new  crystals 

Nj  total  secondary  nucleation 

Nil  Nusselt  number 

Niij  turbulent  Nusselt  number 

Nu0  Nusselt  number  for  a  particle  in  a  station¬ 

ary  fluid 

Pe  Peclet  number 

Pr  Prandtl  number 

qrb  frequency  of  collisions  between  crystals 

and  boundary 

q(r i,r2)  frequency  of  collisions  between  crystals 

of  size  r\  and  r2 
Q*  net  heat  transfer 

r  major  linear  dimension  of  frazil  crystals 

q  face  dimension  of  disk 

R  region  of  phase  space 

Rc  collision  radius 

Re  Reynolds  number 

Rh  hydraulic  radius 

SN  (h)(h) 


t  time 

T f  temperature  of  fluid 

u'  turbulent  fluctuation  of  velocity 

Urb  relative  velocity  of  crystal  and  boundary 
v(ri,r2 )  relative  velocity  of  crystals  of  size  q 

and  r2 

Ve  ( R,t )  external  phase  space  convective  velocity 

Vjf  fluid  velocity 

V;  (R,t)  internal  phase  space  convective  velocity 

V  particle  phase  space  velocity 

V  ( R,t )  phase  space  velocity 

a  thermal  diffusivity 

aT  turbulent  intensity 

8 (r-rc)  Dirac  delta  function 

collision  efficiency 
e  turbulent  dissipation  rate 

q  dissipation  length  scale 

0  bulk  supercooling 

v  fluid  kinematic  viscosity 

pf  density  of  water 

p  £  mass  concentration  of  water  in  mixture 

Pi  mass  concentration  of  ice  in  mixture 

Pi  density  of  ice 

£m  length  scale  of  maximum  eddy  size 

INTRODUCTION 

In  this  chapter,  we  will  discuss  the  development 
of  a  physically  based  quantitative  model  that  de¬ 
scribes  the  dynamic  evolution  of  the  frazil  crystal 
size  distribution  function  during  the  first  stage  of 
frazil  evolution  in  natural  water  bodies.  Two  equa¬ 
tions  serve  as  the  basis  for  the  model:  the  crystal 
number  continuity  equation  and  the  heat  balance 
for  a  differential  volume.  The  use  of  only  two  equa¬ 
tions  simplifies  the  presentation  but  is  not  meant  to 
suggest  any  limitation;  conservation  of  impurities 
(such  as  salt),  for  example,  could  be  added  (see 
Chapter  5).  The  crystal  growth  rate  and  secondary 
nucleation  rate  are  the  major  parameters  that  appear 
in  these  equations.  Expressions  for  both  are  derived. 
We  will  see  that  there  are  three  basic  environmental 
parameters:  1)  the  heat  loss  rate  and  2)  the  crystal 
seeding  rate,  which  appear  directly  in  the  basic 
equations,  and  3)  the  turbulent  energy  dissipation 
rate,  which  influences  both  the  crystal  growth  rate 
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and  the  rate  of  secondary  nucleation.  One  material 
property,  the  number  of  crystals  produced  per  unit 
of  collision  energy,  directly  controls  the  rate  of  sec¬ 
ondary  nucleation. 


BASIC  EQUATIONS 

Crystal  number  continuity  equation 

The  crystal  distribution  will  be  described  in  a 
space  termed  the  crystal  phase  space  or,  more  gener¬ 
ally,  the  particle  phase  space.  Particle  phase  space  is 
defined  by  the  least  number  of  independent  coordi¬ 
nates  that  provides  a  complete  and  useful  descrip¬ 
tion  of  the  properties  of  the  crystal  distribution.  It  is 
convenient,  if  somewhat  arbitrary,  to  divide  particle 
phase  space  into  two  subregions  defined  by  external 
coordinates  and  internal  coordinates.  The  external 
coordinates  describe  the  spatial  distribution  of  the 
crystals.  Internal  coordinates  refer  to  properties  at¬ 
tached  to  each  individual  crystal,  which  quantita¬ 
tively  measure  its  state,  and  are  independent  of  its 
position. 

To  begin,  we  will  consider  a  crystal  distribution 
function  n(R,t).  This  function  is  defined  over  a  re¬ 
gion  R  of  the  particle  phase  space  consisting  of  the 
three  spatial  dimensions  (the  external  coordinates) 
plus  any  number  of  internal  property  coordinates. 
In  all  further  cases  the  internal  coordinates  will  be 
restricted  to  one,  which  will  correspond  to  a  major 
linear  dimension  r  of  the  ice  crystals.  The  function 
n(R,t)  is  defined  as  the  population  density  of  crystals 
in  the  region  R.  At  a  time  t  the  number  dN  of  crystals 
in  an  incremental  region  dR  of  the  particle  phase 
space  is  given  by 

dN  =  ndR  (H) 

and  the  total  number  of  crystals  in  a  region  R  at  time 
t  is 

N(t)  =  l  ndR-  (12) 

R 

It  may  be  necessary  to  deal  with  the  sudden  appear¬ 
ance  (birth)  or  disappearance  (death)  of  crystals  per 
unit  region  in  the  particle  phase  space.  The  net  ap¬ 
pearance  in  an  incremental  region  dR  at  a  time 
would  be 


(B-D)dR 


The  population  balance  of  crystals  in  some  fixed 
region  R,  which  moves  convectively  with  the  par¬ 
ticle  phase  space  velocity  V,  can  be  defined  as 


j-JRndR  =  j(B-D)dR  .  (14) 

Expanding  the  first  term  using  Leibnitz's  rule,  not¬ 
ing  that  the  region  R  was  arbitrary,  and  that  only 
a  single  internal  coordinate  r  is  considered,  where 
G(r,t)  is  the  convective  velocity  along  r  or  simply  the 
growth  rate  of  the  ice  crystal,  we  then  see  that 

—  + — (Gn)  +  D-B  +  V  (ven)  =  0  (15) 

where  Ve  is  the  external  phase  space  convective  ve¬ 
locity. 

This  is  the  number  continuity  equation  in  general 
form.  Further  extensions  can  be  made  by  consider¬ 
ing  the  diffusion  of  crystals,  and  the  rise  velocity  of 
the  crystals. 


Heat  balance 

The  general  expression  for  the  heat  balance  of  the 
frazil-ice-water  system  will  be  developed  strictly  for 
frazil  crystals  suspended  in  fresh  water. 

Consider  a  differential  volume  in  which  pf  is  the 


mass  concentration  of  water  (grams  of  water  per  cu¬ 
bic  centimeter  of  mixture),  p,  is  the  mass  concen¬ 
tration  of  ice  and  the  temperature  of  the  water  is  Tf. 
It  is  assumed  that,  to  good  approximation 


Pf  +  Pi  =  Pf  s  pf 
A  second  assumption  is  then 
Qpi/Cpf  «1 

where  C  is  the  heat  capacity  of  water  and  Q  is  the 
heat  capacity  of  ice. 

Additionally,  heat  conduction  can  be  neglected 
and  heat  capacities  and  the  latent  heat  can  be  con¬ 
sidered  constant  because  of  the  small  variation  in 
temperature.  Therefore,  based  on  the  above  as¬ 
sumptions,  the  heat  balance  for  a  mixture  of  frazil 
and  water  can  be  written 


where  B(R,t)  and  D(R,t)  represent  birth  and  death 
functions  at  a  point  in  the  phase  space. 
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— \—  Pi  +  v(vePi)l  +  7r- 
CpfLdf  V  eKl'J  Cpf 


(16) 


where  Vf  =  convective  velocity  of  the  fluid 

Q*  =  net  heat  transfer  from  the  mixture 
volume 

L  =  latent  heat  of  fusion  of  the  ice  at  the 
equilibrium  temperature  of  the  mix¬ 
ture. 

Parameters  in  the  basic  equations 

The  two  basic  equations  are  those  for  the  crystal 
number  continuity  (eq  15),  and  the  heat  balance  (eq 
16).  The  various  parameters  that  appear  in  these  two 
equations  will  now  be  discussed. 

1.  G — Along  the  a-axis  of  frazil  ice  disk,  it  is  as¬ 
sumed  that  G  is  effectively  determined  by  the  heat 
transfer  rate.  Thus,  in  general 

G  =  ^e 
PiL 

where  h,  the  heat  transfer  coefficient,  is  a  function  of 
the  crystal  size  r  and  the  level  of  turbulence  e,  and  9 
is  the  supercooling  of  the  mixture. 

2.  D — The  death  function  can  be  set  to  zero  for  all 
sizes  of  crystals.  This  is  equivalent  to  assuming  that 
there  is  no  large-scale  breakage  of  the  crystals. 

3.  B — The  birth  function  is  determined  by  the 
rate  of  the  sudden  appearance  of  new  crystals.  New 
crystals  can  appear  as  a  result  of  spontaneous  nucle- 
ation,  secondary  nucleation  and  the  introduction  of 
crystals;  however,  spontaneous  nucleation  is  not 
possible  under  frazil-forming  conditions.  Therefore, 
B  will  be  determined  by  the  rate  at  which  new  crys¬ 
tals  are  introduced  and  the  rate  of  secondary  nucle¬ 
ation. 

Let  Nj  be  the  rate  of  secondary  nucleation.  Nj 
is  a  function  of  the  crystal  size  distribution  n,  the  tur¬ 
bulence  dissipation  rate  e,  the  supercooling  of  the 
mixture  0  and  perhaps  other  parameters.  Let  Nj  be 
the  rate  at  which  new  crystals  are  introduced.  We 
assume  that  new  crystals  are  created  and  introduced 
at  a  size  equal  to  the  critical  radius  of  the  crystals  rc. 
Thus 

B  =  [Nj  (9,  m,  e)  +  Nj  ]  8(r  -  rc ) 


5.  Q* — This  is  determined  by  the  environment  of 
the  water  body  of  interest  and,  in  particular,  the  me¬ 
teorological  and  hydraulic  conditions. 

6.  Ve,  Vf  — The  convective  velocity  of  the  ice  crys¬ 
tals  and  the  fluid  will  generally  be  very  similar.  The 
action  of  buoyancy  and  inertial  forces  on  the  ice 
crystals  may  cause  the  ice  crystal  velocity  to  differ 
from  the  fluid  velocity  if  these  forces  become  large 
compared  to  the  fluid  drag  force. 

Substituting  the  above  expressions  into  eq  15  and 
16  gives  us 


dn  0  9 
9f  p;L  dr 


(hn)  +  v(venj 


=  (NT  +  Nl)8(r-rc) 


(17a) 


Ipi^v 

Cpf 


+  V  Vj  r^ndr  +  — — 
o  *-Pf 


where  n  =  n(x,y,z,r,t) 

0  =  B(x,y,z,t) 
h  =  h(r,e) 

Ve  =  Ve(x,y,zrr,t) 

Nj  =  Nj(d,£,n) 

Vf  =  Vf(x,y,z,t) 

Q*  =  Q*(x,y,z,t). 

Writing  the  equations  in  this  form  emphasizes 
the  dynamic  way  in  which  they  interact.  To  deter¬ 
mine  r  and  n  uniquely,  both  equations  must  be 
solved  simultaneously,  and  the  boundary  condi¬ 
tions  and  initial  conditions  of  r  and  n  must  be 
known.  Difficulties  arise  because  0  and  n  are  di¬ 
mensionally  incompatible. 


where  8 (r—rc)  is  the  Dirac  delta  function. 

4.  pj  — The  mass  of  ice  per  unit  volume  of  mix¬ 
ture  can  conveniently  be  determined  using  the  mo¬ 
ment  equation 

Pi=PiM  n(r)dr 

where  Kvis  a  crystal  shape  function  (Kv  =  An/ 3  for 
the  sphere,  1  for  a  cube). 


ICE  CRYSTAL  GROWTH  RATES 

Intrinsic  kinetic  growth  rate 

The  mechanisms  that  determine  the  rate  at  which 
an  ice  crystal  can  grow  are  the  transport  of  water 
molecules  to  the  crystal  surface  (for  ice  grown  in 
pure  water  this  is  not  a  consideration),  their  incorpo¬ 
ration  into  the  crystal  surface  and  the  transport  of 
latent  heat  away  from  the  surface.  The  incorporation 
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of  molecules  is  controlled  by  the  crystallization  or 
interface  kinetics  of  ice,  and  the  heat  transfer  reflects 
the  particular  physical  situation  of  the  system  under 
consideration. 

The  interface  kinetics  of  ice  has  been  studied  both 
theoretically  and  experimentally.  As  noted  in  Chap¬ 
ter  1,  ice  has  two  principal  growth  directions:  fl-axis 
and  c-axis.  The  interface  kinetics  of  each  growth  di¬ 
rection  appears  to  be  different.  Growth  in  the  c-axis 
probably  proceeds  by  surface  nucleation  for  perfect 
crystals  and  by  a  dislocation  mechanism  for  dam¬ 
aged  crystals.  The  interface  kinetics  of  c-axis  growth 
has  not  been  completely  defined;  the  mechanism  is 
probably  that  of  continuous  growth.  However,  it  ap¬ 
pears  that  the  kinetics  is  very  fast  and  that  for  practi¬ 
cal  purposes  the  growth  rate  of  the  c-axis  is  totally 
controlled  by  the  rate  at  which  heat  is  transported 
away  from  the  interface.  Growth  along  the  c-axis  is 
much  slower  than  that  along  the  c-axis  for  all  sizes  of 
crystals.  This  implies  that  c-axis  growth  is  controlled 
by  the  intrinsic  kinetics.  The  c-axis  growth  rate  does 
not  appear  in  the  number  continuity  equation.  The 
latent  heat  released  by  growth  along  the  c-axis  may 
contribute  somewhat  to  the  overall  heat  balance; 
however,  it  may  be  that  the  latent  heat  released  by  c- 
axis  growth  is  effectively  negligible.  Therefore,  only 
the  growth  along  the  c-axis  will  be  considered. 

Heat  transfer  from  ice  crystals 
suspended  in  turbulent  water 

In  this  section  expressions  for  the  rate  of  heat 
transfer  from  suspended  ice  crystals  will  be  formu¬ 
lated.  To  determine  the  transfer  rate,  it  is  necessary 
to  describe  the  ambient  velocity  distributions  of  the 
fluid  about  the  crystal.  Frazil  is  created  and  devel¬ 
ops  only  in  water  that  is  turbulent.  Rivers  and  chan¬ 
nels  are  inherently  turbulent  because  of  the  instabil¬ 
ity  of  their  bulk  currents,  and  wind  can  make  large 
water  bodies  turbulent.  Frazil  is  also  created  in 
crystallizers  in  which  the  water  is  made  turbulent 
by  impellers,  turbines  or  other  means.  To  describe 
the  velocity  distribution  of  the  water  surrounding 
the  crystals  requires  knowledge  of  the  properties 
and  characteristics  of  turbulence. 

Turbulence  can  be  visualized  as  numerous  inter¬ 
acting  eddies  of  all  possible  scales.  The  very  largest 
eddies  originate  directly  from  the  instabilities  of  the 
mean  bulk  flow.  The  scale  and  orientation  of  these 
largest  eddies  are  imposed  by  the  geometry  of  the 
flow  situation.  Energy  is  extracted  from  the  large  ed¬ 
dies  through  the  inertial  interaction  of  these  eddies 
with  smaller  eddies.  The  energy  cascade  is  not  af¬ 
fected  by  the  fluid  viscosity  until  the  smallest  scales 
are  reached,  where  this  energy  is  dissipated  by  the 


viscosity.  The  dissipation  rate  must  equal  the  rate  at 
which  energy  is  supplied  to  the  small-scale  eddies. 
The  dissipation  rate  8  and  the  fluid  dynamic  viscos¬ 
ity  v  form  a  length  scale  such  that 

r|~(v3/e)1/4  (18) 

where  r|  is  the  dissipation  length  scale  or  the  Kolmo¬ 
gorov  scale. 

If  the  crystal  size  is  small  relative  to  the  Kolmo¬ 
gorov  length  scale,  it  is  in  the  dissipative  regime.  In 
the  dissipative  regime  the  fluid  eddies  are  strongly 
dampened  and  dissipated  by  the  fluid  viscosity.  In 
effect,  the  crystal  is  smaller  than  the  smallest  scales 
of  the  turbulent  eddies.  It  does  not  experience  the 
turbulence  as  interacting  eddies  but  rather  as  a  fluid 
motion  that  varies  linearly  with  position. 

If  the  crystal  size  is  large  relative  to  the  Kolmo¬ 
gorov  length  scale,  it  is  in  the  inertial  regime.  The 
ambient  velocity  can  be  characterized  in  many 
different  ways,  each  corresponding  to  a  different 
eddy  size.  It  seems  reasonable  to  assume,  follow¬ 
ing  Wadia  (1974),  that  the  predominant  shear  that 
the  particle  will  experience  will  be  produced  by 
eddies  closest  to  the  particle  that  are  of  the  same 
size  as  the  particle.  Eddies  that  are  significantly 
larger  than  the  particle  will  entrain  both  the  par¬ 
ticle  and  the  fluid  around  it.  Very  small  eddies 
relative  to  the  particle  size  may  enhance  the  over¬ 
all  transport  by  some  mechanism  of  renewal  of 
the  boundary  layer  surrounding  the  particles,  but 
eddies  of  a  size  comparable  to  the  size  of  the  par¬ 
ticle  will  cause  the  most  significant  gradients  near 
the  crystal  surface. 

The  small-scale  motion,  smaller  than  the  size  of 
the  crystal,  may  enhance  the  heat  transport  from 
the  crystal  by  penetrating  the  boundary  layer 
around  the  crystal.  It  is  difficult  to  quantify  this 
process  but  this  enhancement  has  been  success¬ 
fully  (although  empirically)  accounted  for  by  cor¬ 
relation  with  the  turbulent  intensity  aT  of  the 
fluid,  aj  is  defined  as 


ct-T  =  "^Mf2 /V  (19) 

where  u '2  is  the  rms  valueof  the  velocity  devia¬ 
tion  from  the  mean  velocity  Vf . 

A  method  of  determining  the  Nusselt  number 
that  provides  an  intuitively  easier  means  of  seeing 
the  relative  value  of  the  actual  heat  transfer  coeffi¬ 
cient  is  as  follows.  Let  Niij  be  the  turbulent  Nusselt 
number  defined  as 
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Figure  6.  Nusselt  number  relationships. 
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Let  m*  =  r/r\.  The  heat  transfer  relationships  are  then 
as  follows: 

For  m*<  ly(Pr)1/2 

Nut  =  (1/m*)  +  0.17  Pr1'2  ■  (20a) 

Forl/(Pr)1/2<m*<  =  10 

Nuj  =  [(1/m*)  +  0.55  (Pr/m*)^3] .  (20b) 

For  m*  >  1,  with  a  low  intensity  aTm*4/3  <  1000 


NuT=[(l/m*)  +  0.70a$035(Pr/m)1/3] .  (20c) 

For  m*  >  1,  with  a  high  intensity  a Tm*4/3  >  1000 


Nmt  =  (1/  m* )  +  0.70  a  j25  Pr1/3  (20d) 

These  Nusselt  number  relationships  are  shown  in 
Figure  6. 

Secondary  nucleation 

The  processes  that  govern  the  rates  of  secondary 
nucleation  are  poorly  understood.  However,  a  par¬ 
tial  modeling  of  the  kinetics  of  secondary  nucleation 
is  possible  based  on  the  work  of  Evans  et  al.  (1974a, 
b),  who  demonstrated  that  for  ice  the  production 


rate  of  potential  nuclei  of  new  crystals  and  their  re¬ 
moval  from  the  parent  crystals  could  be  uncoupled. 
The  most  widely  accepted  source  of  the  potential 
nuclei  is  surface  irregularities  that  are  sheared  from 
the  surface  of  the  parent  crystals  (microattrition). 
Two  general  mechanisms  of  removal  of  the  nuclei 
from  the  surface  of  the  parent  crystals  have  been 
suggested:  collisions  of  the  crystals  with  hard  sur¬ 
faces  (including  other  crystals)  and  fluid  shear.  If  the 
rate  of  secondary  nucleation  is  limited  by  the  pro¬ 
duction  rate  of  potential  nuclei,  increases  in  the 
number  of  collisions  of  an  individual  crystal  will  not 
increase  the  production  of  new  crystals.  If  the  rate  of 
secondary  nucleation  is  removal-limited,  however, 
the  parent  crystals  will  produce  the  same  number  of 
new  nuclei  each  collision,  independent  of  the  crys¬ 
tal's  time  history.  From  their  experimental  work, 
Evans  et  al.  (1974a,b)  concluded  that  the  secondary 
nucleation  of  ice  was  limited  by  the  rate  at  which 
potential  nuclei  were  removed  from  the  crystal  sur¬ 
face.  Therefore,  it  was  possible  to  determine  the 
overall  nucleation  rateNj,  with  two  or  more  mech¬ 
anisms  of  removal,  as  the  linear  sum  of  the  actual 
nucleation  rate  attributable  to  each  mechanism  of 
removal  (Nj) 

NT=N1  +  N2  +  ...Ni  (21) 

The  nucleation  rate  of  each  mechanism  of  collision 
can  be  expressed  as  the  product  of  three  functions 
(Botsaris  1976) 

NT=(Et)(F1)(P2)  (22) 
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where  Et  =  rate  of  energy  transfer  to  crystals  by 
collision 

Fj  =  number  of  particles  generated  per  unit 
of  collision  energy 

F2  =  fraction  of  particles  surviving  to  be¬ 
come  nuclei. 

At  this  time  the  values  of  Fj  and  F2  must  be  deter¬ 
mined  empirically.  Therefore,  to  simplify  matters, 
let  Fx  and  F2  be  combined  and  eq  22  be  rewritten  as 

Nt  =  EtSN  (23) 

where  SN  =  (Fi)(F2). 

We  expect  drat  SN  is  a  function  of  all  the  param¬ 
eters  that  govern  the  surface  morphology  and  the 
crystal  growth,  including  supercooling  q,  impurity 
concentrations  Op,  turbulence  level  e,  etc.  The  total 
nucleation  rate  can  be  expressed  as 

Nt  =SN(0,e,Cx,  etc.)  (Etl  +  Et2  +  Et3 . . . j .  (24) 

SUMMARY 

In  this  chapter  a  physically  based  quantitative 
model  of  frazil  ice  was  described.  Two  equations 
serve  as  the  basis  of  the  model:  the  crystal  number 
continuity  equation  and  the  heat  balance  for  a  differ¬ 


ential  volume.  One  focus  of  the  chapter  has  been  de¬ 
scribing  in  detail  two  of  the  major  parameters  that 
appear  in  this  model:  the  crystal  growth  rate  and  the 
secondary  nucleation  rate.  Both  involve  complex  in¬ 
teraction  of  the  crystals  with  the  fluid  turbulence. 
Three  basic  environmental  parameters  and  one  ma¬ 
terial  property  control  the  evolution  of  the  crystal 
size  distribution  function.  These  environmental  pa¬ 
rameters  are  the  heat  loss  rate,  the  turbulent  energy 
dissipation  rate  of  the  fluid,  and  the  crystal  seeding 
rate;  the  material  property  is  the  number  of  crystals 
produced  per  unit  of  collision  energy.  A  tremendous 
amount  of  literature  describes  the  first  two  of  these 
parameters,  and  virtually  none  the  third;  a  modest 
literature  is  available  on  secondary  nucleation,  with 
very  little  describing  ice. 

The  importance  of  a  quantitative  model  such  as 
described  here  is  that  it  allows  predictions  to  be 
made  that  can  be  tested  in  laboratory  experiments. 
By  controlling  or  measuring  the  environmental  pa¬ 
rameters  in  experiments,  knowledge  of  the  resulting 
size  distributions  would  allow  estimates  of  the 
growth  rate  and  secondary  nucleation  rate  to  be 
tested.  Unfortunately,  the  measurement  of  the  crys¬ 
tal  size  distribution  is  quite  difficult,  even  in  the 
laboratory.  However,  such  tests  would  advance  our 
fundamental  understanding  of  frazil  ice  evolution, 
and  the  efficient  solution  of  practical  problems 
could  not  be  far  behind. 
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CHAPTER  5 

Numerical  Simulation  of  Frazil  Ice 

ANDERS  OMSTEDT 


NOMENCLATURE 


Cf 

Cg 

C^Cie,C2£,C3£ 

de 

d{ 

ds 

E 

m 

MO 

f 

GQ>GColl>Gs.GH 

k 

kw 

L 

Nu 

P 

Pb 

Ps 


q 
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T 


Tf 


T 

1n 

t 

u 


«r 

V_ 

V 


x 


y 


frazil  ice  concentration,  volume 
fraction 

grease  ice  concentration,  volume 
fraction 

constants  in  the  turbulence  model 
specific  heat  of  sea  water 
diameter 
crystal  thickness 
crystal  face  diameter 
diameter  for  a  sphere 
coefficient  for  collision  efficiency 
mass  exchange  or  seeding 
net  heat  loss  at  the  air/ water  in¬ 
terface 

Coriolis'  parameter 
source-sink  terms  attributable  to 
ice  formation 
turbulent  kinetic  energy 
thermal  conductivity  of  water 
latent  heat  of  ice 
Nusselt  number 
pressure 

source-sink  term  attributable  to 
buoyancy  production 
source-sink  term  attributable  to 
shear  production 
heat  transfer  between  ice  and 
water 

water  salinity 
water  temperature 
freezing  temperature 
temperature  of  maximum  den¬ 
sity 

maximum  supercooling 
supercooling  at  seeding 
time 

characteristic  time 

mean  velocity  in  x-direction 

collision  velocity 

mean  velocity  in  y-direction 

mean  velocity  vector 

horizontal  coordinate,  positive 

in  east  direction 

horizontal  coordinate,  positive 
in  north  direction 


z  vertical  coordinate,  positive  upwards 

e  dissipation  of  turbulent  kinetic  energy 

v  kinematic  viscosity 

Vj  kinematic  turbulent  viscosity 

pf  frazil  ice  density 

Pg  grease  ice  density 

pM  mixture  density 

pw  water  density 

Oq  turbulent  Schmidt  number  for  Cf 

oc  turbulent  Schmidt  number  for  C„ 

Ok  S  turbulent  Prandtl-Schmidt  number  for  k 

ag  turbulent  Schmidt  number  for  S 

an  turbulent  Prandtl  number  for  heat 

ae  turbulent  Prandtl-Schmidt  number  for  e 


INTRODUCTION 
Scope  of  the  chapter 

The  purpose  of  this  chapter  is  to  review  the  nu¬ 
merical  modeling  of  frazil  ice  in  water  bodies.  The 
recent  developments  in  the  field  of  computational 
fluid  mechanics  have  made  it  possible  to  include 
advanced  turbulence  models  in  environmental  flow 
simulations.  This,  together  with  increased  knowl¬ 
edge  about  the  thermal  regime  during  freezing  and 
frazil  ice  dynamics,  provides  us  with  a  good  starting 
point  for  numerical  simulations  of  frazil  ice. 

To  begin  with  we  should,  however,  recall  that 
models  do  not  give  any  new  information.  Instead 
they  provide  us  with  a  logical  framework  in  which 
we  can  add  our  knowledge  and  test  our  ideas — as 
such  a  "laboratory,"  in  which  we  do  not  get  cold  and 
wet. 

Physical  processes  to  be  considered 

Initially,  ice  forms  in  turbulent  waters  by  the 
introduction  of  ice  crystals  from  the  atmosphere, 
which  become  suspended  in  the  supercooled  water 
(see  Fig.  7).  As  ice  formation  proceeds,  the  sus¬ 
pended  ice  crystals  increase  in  size  and  number, 
stick  to  the  bottom  or  rise  towards  the  surface,  floc¬ 
culate  and  form  clusters  with  high  porosity.  In  the 
ocean,  salt  is  also  rejected  from  the  ice  crystals.  The 
physical  processes  to  be  considered  therefore  in¬ 
clude  cooling,  turbulent  mixing,  mass  exchange 
with  the  air  and  the  water,  and  ice  growth,  multipli- 


25 


Supercooling,  seeding  Frazil  ice  formation,  growth,  Flocculation,surface  ice  and 

break-up,  multiplication,  salt-  anchor  ice  formation 
rejection 


Figure  7,  Process  of  initial  ice  formation  in  turbulent  water  bodies. 


cation  and  flocculation.  The  frazil  ice  evolution  is 
thus  a  complex  interaction  among  fluid  flow,  heat 
and  mass  transfer. 


GUIDANCE  AND 
SIMPLIFYING  ASSUMPTIONS 

Turbulence 

Frazil  ice  forming  under  turbulent  conditions  im¬ 
plies  that  we  need  to  consider  a  relevant  boundary 
layer  flow.  In  oceans,  a  turbulent  Ekman  flow,  and  in 
rivers,  a  turbulent  channel  flow  are  the  most  basic 
boundary  layers.  These  boundary  layers  will  there¬ 
fore  be  discussed  here. 

Different  approaches  to  solving  the  boundary 
layer  equations  are  available,  yet  they  all  start 
from  conservation  equations  for  mean  properties, 
e.g.,  mean  temperature  and  mean  momentum. 
The  different  approaches  adopted  can  basically  be 
divided  into  integral  models  and  closure  models 
of  turbulence.  The  integral  or  mixed  layer  models 
are  based  upon  the  observation  that  the  mean 
properties  are  often  time-dependent  but  constant 
in  a  vertical  direction  within  a  mixed  layer.  Verti¬ 


cal  integration  of  the  conservation  equations  then 
gives  ordinary  differential  equations  for  the  bulk 
properties.  The  closure  models  of  turbulence  in¬ 
volve  the  Reynolds  averaging  of  the  conservation 
equations  and  the  derivation  of  the  average  prod¬ 
ucts  from  different  assumptions  about  the  mixing 
processes. 

The  approach  adopted  in  this  chapter  is  to 
choose  a  closure  model  of  turbulence  with  the  clo¬ 
sure  level  according  to  a  two-equation  model  of  tur¬ 
bulence.  This  kind  of  model  has  been  tested  with 
success  for  many  different  kinds  of  flows.  The 
model  therefore  owes  its  strength  to  its  generality 
(Rodi  1980, 1987). 

Mass  exchange 
and  morphology 

Several  field  observations  report  that  the  frazil 
ice  formation  begins  at  supercoolings  on  the  order  of 
0.01-0.1  °C.  The  initial  nucleation  of  frazil  ice  in  na¬ 
ture  is  therefore  believed  to  be  caused  by  seeding  or 
mass  exchange  with  the  atmosphere  (Osterkamp 
1977).  In  laboratory  studies,  we  can  also  observe 
that,  for  temperatures  close  to  freezing,  the  amount 
of  supercooling  before  frazil  ice  starts  to  form  is  al- 
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most  completely  determined  by  the  seeding  time 
(Hanley  and  Tsang  1984).  Field  and  laboratory  ob¬ 
servations  also  indicate  that  the  initial  crystals  grow 
into  small  disks  and  that  further  ice  production  is 
associated  with  crystal  multiplication.  These  obser¬ 
vations  guide  us  to  include  mass  exchange  with  the 
atmosphere  and  to  regard  the  ice  as  small  crystals 
suspended  in  the  water. 

Gravity 

Another  obvious  observation  is  that  the  ice  crys¬ 
tals  are  lighter  than  the  water.  Turbulent  kinetic  en¬ 
ergy  is  therefore  needed  to  suspend  the  crystals.  The 
conversion  of  turbulent  kinetic  energy  into  potential 
energy  has  thus  to  be  considered  in  the  modeling 
work.  This  is  a  key  feature  for  the  mixing  of  the  up¬ 
per  layers  of  the  ocean,  and  also  when  studying 
suspended  sediment  in  a  bottom  boundary  layer.  It 
is  also  believed  to  be  important  for  the  frazil  ice  for¬ 
mation. 

Salt  rejection 

During  frazil  ice  formation,  the  temperatures  are 
close  to  the  temperature  of  freezing,  and  no  salt  can 
therefore  crystallize.  Instead,  all  salt  is  probably  re¬ 
jected  from  the  interior  of  the  crystal  out  to  the  crys¬ 
tal  surface,  where  it  is  mixed  with  surrounding  wa¬ 
ter.  In  fact,  freezing  can  work  as  an  effective  method 
to  separate  salt  from  water.  This  implies  that  we 
could  assume  that  almost  pure  ice  crystals  are 
formed  during  frazil  ice  formation  and  that  nearly 
all  salinity  is  rejected  out  into  the  surrounding  water 
body. 

Flocculation 

As  cooling  proceeds,  frazil  floes  (grease  ice)  are 
observed  on  the  water  surface.  They  have  then  over¬ 
come  the  turbulent  mixing  through  gravity.  This 
observation  indicates  that  large  crystals  or  floes  of 
crystals  may  become  so  large  that  mixing  is  insuffi¬ 
cient.  Flocculation  may  therefore  work  as  the  main 
mechanism  for  surface  ice  formation  in  turbulent 
waters. 


REVIEW  OF  THE  LITERATURE 

There  are  a  large  number  of  models  on  mixed 
layer  dynamics  in  the  ocean  and  in  rivers;  however, 
very  few  are  directed  to  the  numerical  simulation  of 
frazil  ice. 

In  the  ocean,  Danard  et  al.  (1983)  presented  a 
model  for  initial  sea  ice  formation.  Frazil  ice  forma¬ 


tion  was,  however,  not  considered  in  their  approach; 
only  cooling  was  considered. 

Bauer  and  Martin  (1983)  illustrated  how  we 
could  treat  ice  formation  in  small  leads.  In  their 
model  they  considered  the  pile-up  of  frazil  floes  to 
result  from  wind  and  waves.  Also  considered  was 
how  the  ice  cover  advanced  upwind  until  the  entire 
lead  was  ice  covered.  Pease  (1987)  examined  ice  pro¬ 
duction  in  a  Bering  Sea  polynya,  and  illustrated  that 
the  polynya  width  is  ascribable  to  a  balance  between 
frazil  ice  production  and  ice  advection.  Bulk  transfer 
coefficients  for  heat  and  momentum  over  leads  and 
polynyas  were  presented  by  Andreas  and  Murphy 
(1986). 

Omstedt  and  Svensson  (1984)  developed  a  frazil 
ice  model  for  a  turbulent  Ekman  flow.  The  basic  idea 
was  that  frazil  ice  should  be  described  by  a  bound¬ 
ary  layer  approach,  in  which  buoyancy  effects  from 
frazil  ice  should  be  included.  Omstedt  (1985a)  veri¬ 
fied  the  model  against  laboratory  data  and  extended 
it  to  the  upper  layers  of  the  ocean  (Omstedt  1985b). 

Mellor  et  al.  (1986)  considered  the  turbulent 
boundary  layer  under  drifting  sea  ice  and  the  re¬ 
sponse  to  melting  and  freezing.  In  the  freezing  sim¬ 
ulation,  a  small  amount  of  supercooling  was  mod¬ 
eled  throughout  the  mixed  layer,  and  the  authors 
suggested  that  this  could  be  associated  with  the  for¬ 
mation  of  frazil  ice.  Steele  et  al.  (1989)  introduced  a 
new  scheme  for  the  simulation  of  the  molecular 
sublayers.  In  the  case  of  freezing,  the  amount  of  su¬ 
percooling  was  increased  because  of  the  large  differ¬ 
ence  between  molecular  heat  and  salt  diffusion. 

Freysteinsson  (1970)  discussed  how  frazil  ice  in  a 
river  can  be  calculated  on  the  basis  of  the  thermal 
heat  balance  and  a  regression  analysis  for  the 
amount  of  open  water.  Ginsburg  (1979)  reviewed 
methods  of  forecasting  frazil  ice  in  rivers.  The  re¬ 
view  illustrated  how  conditions  for  frazil  ice  forma¬ 
tion  could  be  presented  in  a  diagram  for  a  specific 
reservoir.  Abramenkov  (1980)  introduced  a  floccula¬ 
tion  equation  for  the  formation  of  frazil  ice  in  rivers. 

Starosolszky  (1981)  reported  on  a  survey  con¬ 
ducted  by  the  World  Meteorological  Organization 
(WMO)  on  frazil  ice  forecasts.  A  system  of  ice  fore¬ 
casting  was  also  outlined. 

Matousek  (1981)  presented  a  mathematical 
model  for  frazil  ice  in  a  river  with  steady  flow.  The 
model  was  based  on  the  thermal  balance  of  a  river 
and  semi-empirical  relations  for  the  ice.  Conditions 
for  frazil  ice  formation  were  formulated  by  Mat¬ 
ousek  (1984)  on  the  basis  of  a  balance  between  rise 
velocity  and  turbulent  mixing  (see  also  Gosink  and 
Osterkamp  1983). 

Daly  (1984)  reviewed  the  frazil  ice  dynamics  and 
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outlined  the  basic  equation  for  the  crystal  number 
continuity. 

Marcotte  and  Robert  (1986)  presented  a  first  at¬ 
tempt  to  model  anchor  ice  in  rivers. 

Nyberg  (1986, 1987)  and  Sahlberg  (1990)  applied 
the  model  of  Omstedt  (1985b)  to  rivers.  In  the  appli¬ 
cations,  frazil  ice  and  frazil-floc  dynamics  were  con¬ 
sidered  as  well  as  the  melting  of  frazil  ice.  Some  of 
the  results  were  earlier  reviewed  by  Omstedt  (1986). 

In  recent  work,  Svensson  and  Omstedt  (1994) 
presented  a  model  of  frazil  ice  dynamics  where  the 
crystal  number  continuity  equation  was  solved  for  a 
well-mixed  jar.  The  following  processes  were  stud¬ 
ied:  initial  seeding,  secondary  nucleation,  gravita¬ 
tional  removal,  growth  ascribable  to  cooling  of  the 
water  volume  and  flocculation-breakup. 

In  the  next  section  the  basic  equations  are  intro¬ 
duced.  In  the  Details  of  Calculations  section,  the  nu¬ 
merical  scheme  used  by  the  present  author  and  his 
colleagues  is  briefly  discussed.  Results  from  the  nu¬ 
merical  simulations  of  frazil  ice  in  the  laboratory, 
ocean  and  river  flows  are  reviewed  in  the  section 
that  follows.  Then,  in  the  Forecasts  section,  that  prob¬ 
lem  is  addressed.  Finally,  a  Discussion  is  given. 


MATHEMATICAL 

FORMULATION 

The  main  part  of  the  water  in  nature  is  ocean  wa¬ 
ter.  In  fact,  fresh  water  constitutes  only  2.7%  of  the 
global  water  budget  (Ashton  1986).  In  the  salinity 
range  from  fresh  water  to  ocean  water,  we  also  have 
brackish  waters,  which  are  often  close  to  the  coast 
and  therefore  in  areas  with  a  high  degree  of  human 
activity.  For  example,  the  Baltic  Sea  and  adjacent 
seas  (the  world's  largest  estuary)  have  salinities 
ranging  from  almost  zero  to  levels  found  typically  in 
the  ocean. 

Salinity,  temperature  and  pressure  influence  sev¬ 
eral  properties  in  the  water.  For  processes  close  to 
the  surface,  pressure  effects  can  often  be  neglected. 
Pressure  effects  will  therefore  not  be  dealt  with  in 
this  chapter.  SI  units  are  used  for  all  variables  except 
temperature,  for  which  degrees  Celsius  are  used. 

Properties  of  water  and  ice 

The  freezing  temperature  Tf  as  well  as  the  tem¬ 
perature  of  maximum  density  Tm  are  functions  of 
salinity  S 

T{  =  -  0.0575  S  +  0.00171  S15  -  0.000215  S2  (22) 
Tm  =  3.982  -  0.2229  S  (23) 


where  S  is  expressed  in  parts-per-thousand.  See  Ap¬ 
pendix  3  in  Gill  (1982),  and  Caldwell  (1978). 

From  the  equations  we  can  see  that  saline  water 
is  freezing  at  lower  temperatures  compared  with 
fresh  water.  In  oceans  the  freezing  temperature  also 
becomes  greater  than  the  temperature  of  maximum 
density,  a  feature  of  great  importance  in  homoge¬ 
neous,  stratified  seas  where  cooling  is  always  associ¬ 
ated  with  convection.  This  is  in  contrast  to  fresh  or 
brackish  water  bodies,  in  which  cooling  is  associ¬ 
ated  with  the  formation  of  a  stable,  stratified  surface 
layer.  The  surface  water  density  (pw)  is  a  function  of 
temperature  and  salinity.  If  nonlinear  effects  are  ne¬ 
glected,  the  equation  reads 

pw  =  999.84  +  0.068  T  +  S  (0.82  -  0.004  T).  (24) 

See  Appendix  3  in  Gill  (1982)  for  the  complete  equa¬ 
tion. 

The  importance  of  temperature  and  salinity  in 
density  calculations  depends  on  the  studied  water 
body.  For  example,  in  estuaries,  which  are  often 
highly  stratified  because  of  a  brackish  surface  layer 
on  top  of  a  more  saline  layer,  the  salinity  almost 
completely  controls  the  stratification. 

Tire  thermal  conductivity  (k^)  is  slightly  depen¬ 
dent  on  temperature  and  salinity.  According  to 
Caldwell  (1974)  the  equation  reads 

ky,  =  0.5719  (1  +  A{F  -  A2T2  -  A4  S)  (25) 

where  Ai  =  0.003,  Tty  =  1.025  x  10-®  and  A4  =  0.00029. 

The  specific  heat  is  a  function  of  salinity  and  tem¬ 
perature.  Again,  disregarding  the  nonlinear  terms, 
the  specific  heat  for  fresh  water  reads  approximately 

cp(S=0)  =  4217.4  -3.7  T .  (26) 

The  corresponding  equation  for  sea  water  reads 

cp  =  cp(S=0)  +  S(-7 .6  +  0.1  T).  (27) 

See  Appendix  3  in  Gill  (1982)  for  the  complete  equa¬ 
tion. 

The  latent  heat  of  fusion  (Lf)  is  a  complex  func¬ 
tion  of  temperature  and  salinity  (see  Yen  1981  or 
Weeks  and  Ackley  1982).  In  the  case  of  frazil  ice  for¬ 
mation,  all  salt  is  probably  rejected  out  from  the 
crystal,  and  pure  ice  crystals  are  formed. 

Heat  equation 

The  starting  point  for  the  heat  budget  is  the  con¬ 
servation  equation  for  heat  energy.  The  derivation 
followed  is  analogous  to  the  derivation  of  the  Rey- 
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nolds  equation;  thus,  heat  is  decomposed  into  a 
mean  and  a  fluctuating  part,  and  the  equation  is 
then  time-averaged.  If  horizontal  mixing  is  not  dealt 
with,  the  equation  reads 

<28» 


obvious  being  a  balance  between  the  ice  formation 
and  the  local  salinity  increase. 

Frazil  ice  equation 

The  derivation  of  the  frazil  ice  equation  follows 
the  same  path  as  above  (see  also  Omstedt  1985b). 
The  equation  reads 


where  z  =  vertical  space  coordinate  positive  up¬ 
wards 

t  =  time  coordinate 
Vj  -  kinematic  eddy  viscosity 
0H  =  Prandtl  number 

Gji  =  source-sink  term  associated  with  freez¬ 
ing-melting. 

The  left-hand  side  denotes  the  total  change  of  heat, 
thus 

~^;PwcpT  =  0^Pwcp^  +  div^VpwCpTj. 

The  boundary  conditions  for  the  heat  equation  have 
to  consider  the  energy  budget  above  the  surface  and 
the  sediment  heat  from  the  bottom.  As  frazil  forms 
in  open  water  areas,  large  amounts  of  ice  can  be  pro¬ 
duced.  This  is  in  contrast  to  columnar  ice  growth, 
where  the  ice  itself  drastically  reduces  the  heat  ex¬ 
change  to  the  atmosphere. 

Before  solving  eq  28,  we  can  notice  some  basic 
balances.  If  the  equation  is  vertically  integrated  and 
advection  is  not  dealt  with,  the  steady-state  solution 
is  a  balance  between  the  net  heat  loss  and  the  release 
of  latent  heat  associated  with  freezing.  If  advection 
becomes  important,  other  basic  balances  may  exist. 
Advection  may,  for  example,  balance  the  net  heat 
loss  and  inhibit  ice  formation. 

Salt  equation 

From  a  conservation  point  of  view,  introducing 
mean  and  fluctuating  salinities  and  time  averaging 
gives 


dCf  d  Vj  dCf 

dt  dz  Gr ,  dz 
V  w 


+  GCf  -  GCon 


(30) 


where  Cf  =  mean  volume  fraction  of  frazil  ice 
0Cf  =  Schmidt  number  for  frazil  ice 
Gq  =  source-sink  term  associated  with 
freezing-melting 

Gcoii  =  sink  term  associated  with  flocculation. 


The  frazil  ice  rise  velocity  is  included  on  the  left  side 
of  eq  30. 

The  boundary  conditions  at  the  air/water  inter¬ 
face  have  to  consider  the  mass  exchange  with  the 
atmosphere.  At  the  water /bottom  interface,  anchor 
ice  may  have  to  be  considered. 

Three  basic  balances  can  be  noticed.  Firstly,  there 
is  a  balance  between  turbulent  mixing  and  rise  ve¬ 
locity  that  controls  whether  surface  ice  may  form  or 
whether  the  crystals  are  kept  in  suspension.  Sec¬ 
ondly,  there  is  a  balance  between  frazil  ice  formation 
and  flocculation.  When  large  amounts  of  crystals  are 
suspended,  the  probability  that  the  crystals  will 
flocculate  increases,  and  thus  only  a  certain  amount 
of  individual  crystals  can  be  kept  in  suspension. 
Thirdly,  there  is  a  balance  between  frazil  ice  forma¬ 
tion  and  anchor  ice  production.  Almost  all  crystals 
may  in  some  situations  form  ice  on  the  bottom,  and 
no  surface  ice  will  then  form. 


Frazil-floc  equation 

The  approach  adopted  by  the  present  author  is  to 
introduce  a  second  ice  equation  to  simulate  the 
frazil-flocs  (grease  ice).  The  derivation  follows  the 
same  ideas  as  above 


ds  =_a_rvx  as"i+G 

dt  3z^os  dz )  s 


(29) 


dCg 

dt 


_3_ 

dz 


+  GColl 


(31) 


where  Gs  is  the  Schmidt  number  for  salinity  and  Gs 
the  source-sink  term  associated  with  freezing-melt¬ 
ing. 

The  boundary  conditions  have  to  consider 
evaporation  and  precipitation.  Melting  and  freezing 
enter  the  equation  through  the  source-sink  term. 
Again  we  can  notice  some  basic  balances,  the  most 


where  Cg  =  mean  volume  fraction  of  frazil-flocs 
Gq  =  Schmidt  number  for  frazil-flocs 
Gcoii  -  source-sink  term  associated  with 
flocculation. 

The  frazil-floc  rise  velocity  is  included  on  the  left 
side  of  eq31. 
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Table  2.  Source-sink  terms  associated  with  frazil  ice 
and  frazil-floc  formation. 


Spherical  crystals 


Plate  crystals 


GH[Wm~3] 

Gc/s-JJ 

Gsk-1! 

Gcoiik'1} 

6qCf 

6  qCf 

6  qSCf 

3  EU'C} 

ds 

d,PfL 

4  KL 

4  qCf 

4  qCf 

4  qSCf 

£«rCf 

df 

dfpfL 

d(P„L 

de 

The  boundary  condition  is  that  there  are  no 
fluxes  at  the  vertical  boundaries. 

We  can  see  from  the  equation  that  frazil-floes  are 
assumed  to  form  only  because  of  flocculation  and 
not  through  thermal  processes.  Some  basic  aspects 
of  eq  31  can  be  noticed.  Firstly  because  of  larger  ag¬ 
gregates  of  ice,  the  rise  velocity  will  rapidly  over¬ 
come  the  turbulent  mixing,  and  the  ice  will  float  up 
to  the  surface.  Secondly,  the  formation  of  surface  ice 
may  become  an  outcome  of  the  balance  between  ice 
production  and  advection  (see  also  Bauer  and  Mar¬ 
tin  1983). 

In  the  present  approach,  we  thus  distinguish  be¬ 
tween  two  ice  regimes,  frazil  ice  and  frazil-flocs.  An¬ 
other  way  of  looking  at  the  equations  is  to  say  that 
we  distinguish  between  small  crystals,  which  are 
formed  because  of  thermal  processes,  and  larger  crys¬ 
tals,  which  are  aggregates  of  the  smaller  crystals. 

A  summary  of  the  different  source-sink  terms 
associated  with  ice  formation  is  given  in  Table  2, 
where  the  heat  transfer  q  and  the  collision  velocity  ur 
are  according  to 
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where  Nu  =  Nusselt  number 
Tj  =  ice  temperature 
kw  =  thermal  conductivity  of  water 
v  =  kinematic  viscosity 
e  =  dissipation  of  turbulent  kinetic  energy 
de  =  ice  crystal  thickness 
df  =  ice  crystal  face  diameter. 

For  the  derivation  of  the  source-sink  terms,  see 
Omstedt  (1985b). 


It  should  be  noticed  that  the  source-sink  terms 
associated  with  frazil  ice  formation  are  linearly  de¬ 
pendent  on  the  frazil  ice  concentration,  while  the 
source-sink  terms  associated  with  flocculation  are 
proportional  to  the  square  of  the  frazil  ice  concentra¬ 
tion.  The  underlying  assumption  is  that  the  collision 
velocity  is  independent  of  the  frazil  ice  concentra¬ 
tion. 

Mixture  density  equation 

The  frazil  ice  problem  is  a  typical  two-phase  flow, 
where  ice  is  one  phase  and  water  the  other.  In  the 
field  of  computational  fluid  mechanics,  it  is  today 
possible  to  solve  multiphase  flows.  However,  the 
present  author  believes  that  a  diffusion  or  mixture 
formulation  may  hold  for  many  applications.  From 
the  mixture  formulation,  the  distribution  of  the  dis¬ 
persed  phase  is  calculated  from  scalar  equations, 
such  as  eq  30  and  31. 

To  include  the  concentrations  in  the  mixed  layer 
dynamics,  we  introduce  a  mixture  density  equation 
according  to 

Pm  =  Pw  +  (PrPw)Cf  +  (pg  Pw)^-g  (32) 

where  pf  is  the  frazil  ice  density  and  pg  the  frazil- 
flocs'  density. 

When  a  mixture  density  is  introduced  in  the  dy¬ 
namics,  an  increased  frazil  formation  will  be  associ¬ 
ated  with  the  formation  of  a  stably  stratified  surface 
layer.  The  turbulent  mixing  thus  has  to  work  against 
gravity  and  turbulent  kinetic  energy  is  therefore 
transformed  to  potential  energy. 

Boundary  layer  equations 

Starting  from  the  Navier-Stokes  equation,  intro¬ 
ducing  mean  and  fluctuating  velocities  and  time-av¬ 
eraging  the  equation,  we  derive  the  Reynolds  equa¬ 
tion.  In  this  section  we  shall  reduce  the  equation  and 
only  outline  two  basic  boundary  layers,  namely  the 
turbulent  Ekman  flow  and  the  turbulent  channel 
flow. 

The  Ekman  boundary  layer  equations  read 

f-lhf)*'1'  » 

£-«•.£)-/»  >» 

where  U  and  V  are  mean  velocities  in  horizontal  di¬ 
rection  and/is  Coriolis'  parameter. 
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The  channel  boundary  equation  reads 


du  _  d  f  di n _ i_dp 

dt  dz  V  T  dz  J  pw  dx 


(35) 


where  dP/dx  is  the  horizontal  pressure  gradient. 

The  boundary  condition  for  the  air /water  inter¬ 
face  is  related  to  the  wind  stress,  and  at  the  water/ 
bottom  interface,  the  velocities  are  set  equal  to  zero. 

The  fundamental  balances  are  between  the  verti¬ 
cal  shear,  associated  with  the  wind  stress,  and  the 
earth  rotation  or  the  vertical  shear,  associated  with 
the  bottom  friction,  and  the  horizontal  pressure  dif¬ 
ferences. 


Turbulence  model 

Several  turbulence  models  exist  and,  as  stated 
earlier,  we  have  chosen  to  use  a  two-equation  model 
of  turbulence,  with  one  equation  for  the  turbulent 
kinetic  energy  k  and  another  for  the  dissipation  rate 
of  turbulent  kinetic  energy  e  (see  Svensson  1979  or 
Rodi  1987).  The  turbulent  mixing  processes  are  as¬ 
sumed  to  be  caused  by  current  shear  and  convection 
caused  by  cooling  and  salt  rejection.  The  turbulent 
mixing  processes  are  de-  Xable  3>  Constants  in 
scribed  by  vertical  ex-  the  turbulence  model, 
change  coefficients.  The 

introduction  of  exchange  Constant _ Value 


coefficients  and  gradient 
laws  excludes  counter¬ 
gradient  fluxes.  Break¬ 
ing  waves  and  Lang¬ 
muir  circulation  are  not 
dealt  with.  The  equa¬ 
tions  read 
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de  _  _3_[  vT  dt 
dt  3z^ae  dz 


p  _  vTc?  dPw  ^  Pf  pw  dCf 
b  Po  [  ^  <*Cf  dz 

Pg  —  Pw  dCg 

Op  dz 

Lg  ) 


(39) 


vT  =  C^  (40) 

where  Ps  is  the  production  caused  by  shear  and  Ph 
the  production-destruction  caused  by  buoyancy. 
The  constants  in  the  turbulence  model  are  denoted 
by  Clg,  Cip  C3g  and  C^,  respectively,  and  the  values 
can  be  found  in  Table  3. 


DETAILS  OF 
CALCULATIONS 

All  the  relevant  differential  equations  can  for¬ 
mally  be  written  as 

+  div(V<j>)  =  div(T  grad  <j>)  + 

wherej)  =  dependent  variable 
V  =  water  velocity 
F  =  exchange  coefficient 
S,),  =  source  and  sink  terms. 

The  four  terms  in  this  general  equation  are  the  un¬ 
steady  term,  the  advective  term,  the  diffusion  term 
and  the  source-sink  term.  The  dependent  variable 
can  stand  for  heat,  salinity,  frazil  ice,  frazil-flocs,  mo¬ 
mentum,  turbulent  kinetic  energy  or  the  dissipation 
rate  of  turbulent  kinetic  energy. 

All  conservation  equations  can  thus  be  written  in 
a  similar  way,  which  implies  that  our  equations  can 
be  solved  using  the  same  numerical  method.  The 
boundary  conditions  may  be  specified  in  two  differ¬ 
ent  ways;  either  by  prescribing  the  value  or  by  pre¬ 
scribing  the  flux  of  the  variable  in  question. 

When  discretizing  the  general  differential  equa¬ 
tion,  several  schemes  exist.  The  present  author  and 
his  colleagues  have  used  a  control-volume  formula¬ 
tion.  The  formulation  starts  by  dividing  the  water 
body  into  control  volumes  around  grid  points  and 
integrating  the  equation  over  these  volumes.  The 
discretized  equation  obtained  in  this  manner  thus 
concentrates  on  the  conservation  principle,  and  the 
discretization  makes  physical  interpretations  pos- 


sible.  For  further  discussions  about  the  numerical 
scheme  the  reader  is  recommended  to  consult  the 
book  by  Patankar  (1980). 

RESULTS 

In  this  section  some  of  the  results  by  the  present 
author  and  his  colleagues  are  given.  In  all  calcula¬ 
tions  presented,  two  general  equation  solvers  for 
hydrodynamic  flows  have  been  used,  the  PROBE 
program  by  Svensson  (1979)  for  one-dimensional 
boundary  layer  problems  and  the  PHOENICS  pro¬ 
gram  by  Rosten  et  al.  (1982)  for  the  river  applications. 

Laboratory  simulations 

There  are  today  several  important  laboratory  ex¬ 
periments  on  frazil  ice  formation  available,  e.g., 
those  of  Arakawa  (1954),  Carstens  (1966),  Mueller 
and  Calkins  (1978),  Martin  and  Kauffman  (1981), 
Ettema  et  al.  (1984),  Hanley  and  Tsang  (1984),  Tsang 
and  Hanley  (1985). 

One  main  problem  with  these  experiments  is  that 
the  turbulent  flows,  under  which  the  ice  is  formed, 
are  often  complex  and  difficult  to  simulate.  The  ex¬ 
periments  are  either  conducted  in  a  turbulence  jar 
with  a  vertically  oscillating  grid  or  in  flumes  of  race¬ 
track  shape,  where  the  flow  is  driven  by  a  propeller. 

From  a  modeling  point  of  view,  we  need  experi¬ 
mental  data  from  well-established  boundary  layer 
flows.  The  flume  flow  experiment  probably  meets 
this  requirement  in  the  best  way;  however,  the  tur¬ 
bulence  generated  by  the  propeller,  side  wall  effects 
and  secondary  circulation  complicate  the  flow. 

In  Omstedt  (1985a),  the  laboratory  experiments 


of  Tsang  and  Hanley  (1985)  were  analyzed.  In  the 
numerical  simulation,  the  frazil  ice  crystals  were 
treated  as  thin,  uniform  plates.  The  heat  transfer  be¬ 
tween  the  frazil  ice  and  surrounding  water  was  cal¬ 
culated  from  Fourier's  law  of  heat  conduction,  with 
the  heat  transfer  coefficient  based  upon  a  constant 
Nusselt  number.  The  general  experimental  findings 
were  well  reproduced  by  the  model.  For  example, 
the  time  for  the  temperature  curve  to  reach  90%  re¬ 
covery  from  maximum  supercooling,  as  well  as  the 
maximum  supercooling,  were  numerically  well 
simulated  (Fig.  8). 

The  time  evolution  during  initial  frazil  ice  forma¬ 
tion  was  studied  by  comparing  model  calculations 
with  normalized  data  from  Tsang  and  Hanley 
(1985).  In  Figure  9,  the  time  evolutions  of  the  nor¬ 
malized  frazil  ice  concentrations  during  nine  experi¬ 
ments  are  considered.  The  normalized  parameters 
are  the  dimensionless  time,  defined  as  the  ratio  be¬ 
tween  time  and  characteristic  time,  and  the  dimen¬ 
sionless  frazil  ice  concentration,  defined  as  the  ratio 
between  frazil  ice  concentration  and  characteristic 
concentration.  The  characteristic  time  is  defined  as 
the  time  for  the  temperature  curve  to  reach  90%  re¬ 
covery  from  maximum  supercooling,  and  the  char¬ 
acteristic  concentration  is  defined  as  the  frazil  ice 
concentration  at  the  characteristic  time.  The  calcula¬ 
tions  reproduce  the  data  quite  well,  except  at  the 
very  first  stage  of  the  ice  formation.  This  was  attrib¬ 
uted  to  the  growth  of  the  very  first  small  crystals 
and  was  not  dealt  with  by  the  model. 

Ocean  simulations 

In  the  ocean,  frazil  ice  forms  in  at  least  four  differ¬ 
ent  situations  (Martin  1981): 


Figure  8.  Measured  and  calculated  relationship  between  the  characteristic  time  tc,  the  temperature  of  maximum 
supercooling  ATmin  and  the  amount  of  supercooling  at  seeding  ATn  (from  Omstedt  1985a). 
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Measured  and  calculated  relationship  between  the  normalized  frazil  ice  concentration  Citand  the  normalized  time  t*  (from  Omstedt  1985a), 


ice  thidtress  loti)  ,  Surface  Iweperature  CO  (  &rface  ice  coocenfraficn  (*/-) 


1.  In  open  water  areas,  caused  by  a  net  heat  loss  3.  Below  freezing  sea  ice,  owing  to  the  drainage 

to  the  atmosphere  and  similar  to  the  process  found  of  cold,  dense  brine  associated  with  the  ice  growth, 
in  rivers.  4.  Attached  to  ice  shelves  and  icebergs  and  asso- 


2.  At  the  interface  between  two  fluid  layers,  each 
at  its  freezing  point  and  with  different  salinities. 


ciated  with  rising  sea  water. 

The  importance  of  frazil  ice  for  the  sea  ice  cover 
in  the  Antarctic — the  world's  largest 
frazil  ice  producer — is  further  dis¬ 
cussed  by  Lange  (1990). 

This  chapter  considers  only  frazil  ice 
formation  in  open  water  areas.  The 
general  frazil  ice  equations,  eq  30  and 
31,  with  corresponding  source-sink 
terms  can,  however,  also  be  applied  to 
the  other  situations.  The  pressure  ef¬ 
fects  and  the  relevant  mechanism  of 
nucleation  should  then  also  be  consid¬ 
ered. 

In  the  surface  layers  of  the  ocean,  the 
mixed-layer  dynamics  play  a  very  im- 
1  -Tim'  port  ant  role.  With  weak  winds  and 

*96  (hours)  stratified  or  shallow  waters,  ice  can 


rapidly  form  over  large  areas.  How¬ 
ever,  in  deeper  oceans  ice  may  never 
form  because  of  convection. 

In  Figure  10,  the  effect  of  wind  mix¬ 
ing  on  supercooling  and  ice  formation 
is  illustrated.  If  the  wind  speed  in¬ 
creases,  warm  water  from  deeper  lay¬ 
ers  is  mixed  into  the  surface  water,  and 
the  ice  formation  is  delayed. 

The  effect  of  mass  exchange  on  su¬ 
percooling  and  ice  formation  is  illus¬ 
trated  in  Figure  11.  It  is  interesting  to 
note  that,  compared  with  high  rates  of 
mass  exchange,  low  rates  of  mass  ex¬ 
change  lead  to  higher  supercooling 


and  more  rapid  recovery  times.  This 


Figure  10.  Effect  of  wind  mixing  on  supercooling  and  initial  ice formation 
(from  Omstedt  and  Svensson  1984). 


implies  that  even  if  the  rate  of  mass  ex¬ 
change  is  modeled  too  low  or  too  high, 
the  dynamics  tries  to  "correct"  the 
time-temperature  evolution. 

In  the  model  by  Omstedt  and 
Svensson  (1984),  the  frazil  ice  crystals 
were  treated  as  spheres  with  a  mean 
radius  of  10-3  m,  and  flocculation  was 
not  dealt  with;  instead,  the  calculation 
was  terminated  at  a  high  frazil  ice  con¬ 
centration.  The  model  was  extended 
in  Omstedt  (1985b),  where  the  frazil 
ice  crystals  were  treated  as  thin  plates 
and  flocculation  of  frazil  ice  crystals 
was  considered. 

Two  basic  differences  between  the 
calculations  in  Omstedt  and  Svens- 
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son  (1984)  and  in  Omstedt  (1985b)  were  found. 
Firstly,  the  time-temperature  response  to  frazil  ice 
formation  became  more  rapid  (Fig.  12).  Secondly, 
flocculation  or  grease  ice  formation  started  at  lower 
frazil  ice  concentrations  than  was  assumed  in 


Omstedt  and  Svensson  (1984).  In  fact,  the  calcula¬ 
tions  indicate  that  flocculation  becomes  important 
even  for  frazil  ice  concentrations  below  l%o. 

The  amount  of  open  water,  and  thus  the  area 
where  surface  frazil  ice  can  form,  is  an  intricate  bal¬ 


ance  among  ice  drift,  ice  formation  and 
melting.  Polynyas,  leads  and  the  posi¬ 
tions  of  ice  fronts  are  the  areas  where 
this  balance  can  be  in  favor  of  frazil  ice 
formation.  For  a  review  of  different  as¬ 
pects  of  sea  ice  the  reader  is  referred  to 
Untersteiner  (1986). 

River  simulations 

In  rivers  the  advection  and  the  ge¬ 
ometry  play  an  important  role.  In 
some  river  reaches,  the  velocities  are 
low,  and  shore  ice  can  grow  out  and 
cover  the  reach.  In  other  reaches,  the 
water  flow  is  faster,  and  frazil  ice  may 


form.  Again,  we  have  to  understand 
the  basic  balances  that  form  or  break 
up  the  ice  cover. 

A  frazil  ice  model  for  rivers  was  pre¬ 
sented  by  Nyberg  (1986).  The  simula¬ 
tion  of  frazil  ice  followed  that  by 
Omstedt  (1985b),  with  some  modifica¬ 
tions.  The  modifications  were  that  the 
flow  was  treated  as  a  two-dimensional 
turbulent  channel  flow  and  that  melt¬ 
ing  was  considered  as  well  as  the 
formation  of  anchor  ice.  From  some 
idealized  calculations,  Nyberg  (1986) 
presented  vertical,  as  well  as  down¬ 
stream,  distributions  of  temperature, 
frazil  ice  and  frazil  floes  (Fig.  13). 

The  model  has  been  extended  to  a 
river  reach  in  Sweden  (Fig.  14).  In  this 
application,  the  river  flow  was  treated 
as  a  transient  three-dimensional  chan¬ 
nel  flow,  using  body-fitted  coordi¬ 
nates.  The  number  of  grid  cells  used 
was  1560,  and  the  time  step  extended 
over  a  few  minutes.  A  quasi-steady- 
state  situation  was  reached  after  about 
1  hour  of  simulation,  which  required 
about  20  minutes  of  CPU-time  on  a 
VAX  8600  computer.  Applications  of  the 
model  during  two  winter  periods  are 
further  discussed  by  Sahlberg  (1990). 

The  questions  analyzed  by  the 
model  simulations  were:  When  and 


Figure  11.  Effect  of  mass  exchange  on  supercooling  and  ice  formation  where  will  frazil  ice  form  during  dif- 
(from  Omstedt  and  Svensson  1984).  ferent  environmental  conditions? 
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Figure  12.  Water  cooling  and  initial  ice  formation  in 
an  idealized  ocean  situation  (from  Omstedt  1985a). 
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Figure  13.  Water  cooling  and  initial  ice 
formation  in  an  idealized  river  situation 
(from  Nyberg  1986). 
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Figure  14.  Frazil  ice  formation  in  Stornorrfors — an  example.  The  discharge  is  600  m3/s,  the 
inflow  temperature  is  0.01  °C  and  the  net  heat  loss  is  450  W/m2  (from  Nybergl987).  a — Grid  net, 
b — Calculated  water  temperature  (°C),  c — Calculated  frazil  ice  concentration  (Cj/lCr5). 


FORECASTS 

In  this  section,  the  problem  of  forecasting  frazil 
ice  is  discussed.  Any  forecasting  of  the  onset  of  frazil 
ice  will  rely  on  successfully  predicting  the  cooling 
rate.  In  oceans  and  rivers,  models  for  predicting  wa¬ 
ter  cooling  are  available  (see,  e.g.,  Omstedt  1984, 
1990,  Shen  et  al.  1984).  The  results  from  these  appli¬ 
cations  are  encouraging. 

In  this  discussion  it  is  important  to  distinguish 
between  diagnostic  studies  and  forecast  studies.  In 
general,  the  main  studies  in  the  literature  about 
cooling  are  diagnostic,  meaning  that  the  time-tem¬ 
perature  evolutions  are  simulated  on  the  basis  of 
measured  weather  data.  The  forecast  studies  in¬ 
clude  a  further  complication  associated  with  the  un¬ 
certainty  in  the  weather  parameters. 

For  the  prediction  of  the  onset  of  frazil  ice  for¬ 
mation,  some  diagrams  are  available  (Ginsburg 
1979).  The  diagrams  are  probably  quite  site-spe¬ 
cific  and  should  therefore  be  further  developed. 
One  way  of  doing  this  is  to  use  an  advanced  frazil 
ice  model  and  develop  local  diagrams  for  the  wa¬ 


ter  body  of  interest.  Basic  parameters  in  these  dia¬ 
grams  should  be: 

1.  Upstream  temperature 

2.  Discharge-volume  flow 

3.  Air  temperature 

4.  Wind 

5.  Cloudiness 

6.  Humidity 

7.  Amount  of  open  water. 

The  list  above  indicates  that  somewhat  complicated 
diagrams  have  to  be  constructed,  but  for  specific  lo¬ 
cations  they  may  probably  serve  as  useful  tools. 

The  main  elements  in  a  forecasting  system  are 
the  following: 

1.  A  budget  for  the  heat  exchange  between  the 
air/ water  and  the  water /bottom  interfaces. 

2.  A  hydrodynamic  model,  including  the  tem¬ 
perature  equation. 

3.  A  frazil  ice  model,  including  the  mass  ex¬ 
change  with  the  atmosphere. 

4.  An  ice  cover  model,  predicting  the  amount 
of  open  water. 

The  forecaster  needs  the  following  initial  data: 
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Figure  15.  How  oceanographic  and  meteorological  data  can  be  combined  with 
mathematical  models  (after  Omstedt  1990). 


1.  Initial  temperatures  and  salinities. 

2.  Initial  upstream  temperatures  and  dis¬ 
charge-volume  flow. 

3.  Initial  distribution  of  open  water  areas. 

The  strategy  when  forecasting  the  frazil  ice  for¬ 
mation  will  depend  upon  the  specific  water  body 
studied.  However,  the  forecaster  needs  to  be  well  in¬ 
tegrated  with  the  weather  forecasting  service.  Ow¬ 
ing  to  the  limited  number  of  temperature  and  ice 
cover  data,  etc.,  available,  which  are  necessary  for 
starting  up  any  calculations,  the  forecasts  often  need 
to  be  supplemented  with  updating  calculations. 
These  calculations,  based  on  observed  weather,  wind 
and  ice  conditions,  should  then  produce  the  neces¬ 
sary  starting  values  for  the  forecasts  (Omstedt  1990). 

In  Figure  15,  the  model  structure  of  a  forecasting 
system  for  water  temperatures  applied  to  the  Baltic 
Sea  is  illustrated.  The  model  (PROBE-BALTIC)  is 
used  in  two  modes.  The  first  one,  the  real-time 
mode,  fits  the  oceanographical  and  meteorological 
data  collected  in  the  recent  past  to  the  model.  The 
second  one,  the  forecast  mode,  starts  from  the  real¬ 
time  mode,  and  calculates  the  forecasts. 


DISCUSSION 

The  purpose  of  this  chapter  was  to  review  the 
numerical  modeling  of  frazil  ice  in  water  bodies. 
Few  models  are  available,  but  it  is  believed  that  the 
model  approach  given  by  the  present  author  and  his 
colleagues  forms  a  good  base  for  further  research 
and  also  for  applications  to  different  water  bodies. 
More  efforts  need  to  be  applied  to  modeling  and 
measuring  frazil  ice  and  particularly  to  the  crystal 
number  continuity  equation  (see  further  discussion 
in  Svensson  and  Omstedt  [1994]). 

In  diagnostic  studies,  all  equations  should  be 
used,  and  diagrams  for  specific  areas  could  be  con¬ 
structed.  When  developing  forecast  systems,  the 
choice  between  diagrams  and  numerical  models 
will  be  a  question  of  how  accurate  the  forecast  re¬ 
quirements  are.  The  accuracy  of  any  frazil  ice  fore¬ 
cast  is  probably  as  accurate  as  the  water  cooling 
forecast  allows  it  to  be.  This  may  therefore  guide  the 
choice  of  model  complexity  when  developing  fore¬ 
cast  systems  for  frazil  ice. 
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